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Preface 


The  objective  of  this  investigation  was  to  try  to 


better  understand  the  interactions  of  jet  and  annulus  flows 


in  the  presence  of  a  bluff  body.  These  interactions  are 


very  important  m  combustors  and  injectors,  and  a  greater 


understanding  of  flow  interaction  can  be  ultimately  used  to 


increase  the  efficiency  and  performance  of  these  devices. 


This  investigation,  however,  was  conducted  from  a  scientific 


viewpoint,  not  an  engineering  one. 


A  unique  apparatus  was  used  to  determine  the  nature  of 


interactions .  A 


two-dimensional 


bluff 


body 


surrounded  a  planar  slot  jet  which  exhausted  into  a 


parallel,  moving  stream,  the  annulus  flow.  Flew 


visualisation  provided  qualitative  insight  into  flowfield 


occurrences.  Laser  Doppler  Anenometry  was  used  to  measure 


axial  and  transverse  velocities  for  the  two-dimensional 


configuration.  Light  scattering,  Mia  scattering,  was  used 


to  determine  the  shedding  frequency  of  the  large  vortex 


structures  present  in  the  flow.  The  results  of  these 


experiments  were  compiled,  and  a  description  of  the  jet  and 


annulus  inter action  was  developed.  As  in  any  problem  in 


fluid  mechanics,  there  will  be  some  doubt  concerning  the _ 

ion  For 

applicability  of  the  description  for  conditions  and~RAtI 


geometries  different  than  those  tested.  Further  studies, 


;nced 


then. 


should  test  the  "robustness"  of  the  physical 


catio 


description  provided  herein. 
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ABSIBACT 

A  particular  interaction  of  bluff  body  and  jet  flows  was 
studied  using  a  two-dimensional  slot  jat  (aspect  ratio  -  50) 
surrounded  by  a  two-dimensional  bluff  body.  The  co- flowing 
system  used  air  in  the  annulus  and  nitrogen  in  the  jet. 
Flow  visualization  identified  flow  patterns  containing  large 
vortex  structures.  With  a  fixed  annulus  velocity,  drastic 
changes  in  the  flow's  structure  were  observed  as  jet 
velocities  changed.  Laser  Doppler  Anenoaetry  Measurements 
of  instantaneous  velocities  were  examined  to  characterize 
the  flovfield  details,  and  high  speed  movies  provided 
qualitative  insight  into  flow  dynamics.  The  concept  of  jet 
and  annulus  "fields"  is  introduced.  The  "jet  field"  is  the 
region  above  the  bluff  body  that  extends  six  or  seven  slot 
widths  downstream,  and  the  "annulus  field"  is  the  region 
above  the  jet  field.  The  annulus  flow  dictates  the 
flowfield  dynamics  in  the  annulus  field  up  to  jet-to-annulus 
velocity  ratios  of  1.45.  The  jet  is  responsible  for  the 
formation  of  the  large-scale  structures  and  dominates  the 
flow  in  the  jet  field.  Jet  domination  of  the  entire 
flowfield  is  seen  at  velocity  ratios  of  1.45  and  above.  Mie 
scattering  was  used  to  determine  shedding  frequencies  of  the 
large-scale  structures;  Strouhal  number  based  on  jet 
velocity  and  slot  width  was  a  constant  0.098. 


xii 


Li.  Introduction 


Many  atudias  of  wakes ,  plumes,  and  jets  have  been 
accomplished.  With  the  recent  attention  given  to  coherent 
structures  in  turbulent  flow,  many  studies  have  begun  to 
focus  on  the  deterministic  character  of  the  flow.  It  had 
been  felt  that  turbulent  flows  could  only  be  analyzed  using 
statistical  tools.  Some  people  believe  now,  however,  that 
the  character  of  the  flow  is  completely  determined  by 
occurrences  in  the  near  field,  where  the  coherent  structures 
exist  (List,  1982;  Laufer,  1981;  Hussain,.  1986).  Hussain 
reports  that  coherent  structures  may  be  the  tractable  part 
of  turbulence  and  mcy  contain  most  of  the  essential  physics 
of  turbulence.  Laufer  says  that  the  initial  instability  has 
a  strong  influence  on  subsequent  jet  development.  In  his 
review  article,  List  tells  us  that  there  is  now  overwhelming 
evidence  that  the  initial  growth  of  turbulent  jets  is  a 
direct  result  of  large-scale  motions  generated  at  jet 
boundaries.  This  emphasis  has  only  recently  shifted  to  the 
region  near  the  flow's  point  of  initiation  (near  field). 
For  turbulent  flows,  the  near  field  of  plane  jets  is 
considered  to  extend  from  the  jet  exit  to  z/D  -  30  or  40, 
where  z  is  the  axial  dimension  and  D  the  slit  width.  The 
near  field  of  axisymmetric  jets  is  considered  to  end  at 
approximately  z/D  ■  8,  where  D  is  now  the  jet  exit  diameter. 

Thousands  of  papers  on  turbulence  have  been  written, 
but  most  of  these  efforts  focused  on  the  far  field.  Also, 
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nearly  all  work  has  been  accomplished  for  geometries  that 
are  not  complicated  (i.e.,  jets,  shear  layers,  bluff  bodies, 
etc  are  studied,  but  combinations  like  a  jet  flowing  into  a 
moving  stream  are  relatively  rare) .  Studies  that  did 
investigate  co-flowing  systems  showed  that  the  parallel  flow 
has  a  significant  impact  on  the  jet  flow.  Bradbury  (1965) 
studied  a  turbulent  plane  jet  exhausting  into  a  slow-moving 
parallel  flow  and  found  self-similarity  of  the  flow  was 
achieved  much  further  downstream  when  a  parallel  flow  (even 
a  very  weak  one)  was  used.  Everitt  and  Robins  (1978)  have 
recently  confirmed  this  result.  Their  experiments  on  a 
plane  jet  exhausting  into  a  moving  airstream  proposed  a 
concept  of  "strong"  and  "weak"  jets.  The  weak  jet  is  one 
whose  velocity  is  small  compared  to  the  parallel  flow  and 
results  in  a  slow  development  of  the  flow. 

Some  recent  work  on  various  combustor  configurations 
has  introduced  unique  geometries  which  may  be  used  to  better 
characterize  the  events  in  the  near  and  far  fields  of  jets 
and  co-flowing  systems.  Roquemore  et  al  (1986)  studied  an 
axisymmetric  bluff  body  combustor,  and  Ballal  and  Chen 
(1986)  and  Tankin  (1987)  tested  a  two-dimensional  slot  jet 
surrounded  by  a  bluff  body  (the  configuration  will  be 
referred  to  as  the  2-D  centerbody,  and  it  can  be  seen  in 
Figure  3.1).  This  configuration  represents  the 
two-dimensional  analogue  of  an  afterburner's  fuel  ring 
(Ballal  and  Chen,  1986) .  The  bluff  body  is  used  to  provide 
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stability.  The  high  aspect  ratio  jat  axhauata  into  a  moving 
airstreait,  ra  far  rad  to  as  tha  annulus  flow.  Plane 
stabilisation  is  greatly  enhanced  by  the  fluid  dynamic 
processes  that  occur  in  the  near  field.  This  configuration 
is  indeed  unique  since  it  incorporates  a  jet,  a  bluff  body, 
and  a  parallel  annulus  flow.  Previous  research  on  another 
configuration  like  this  has  not  been  identified.  This 
set-up  is  also  important  in  that  it  provides  a  geometry 
which  simplifies  advanced  Computational  Fluid  Dynamics  (CFD) 
models. 

Belial  and  Chen  studied  a  jet  flow  and  used  C02  as  the 
working  fluid.  Annulus  velocity  was  only  10  percent  of  the 
jat  velocity,  and  this  ensured  they  had  a  jet  flow.  From 
the  Laser  Doppler  Anenometry  (LDA)  data  they  obtained,  which 
focused  on  the  near  field,  they  were  able  to  provide  a 
physical  description  of  the  flowfield  development.  Attempts 
to  verify  this  description  using  flow  visualisation  failed. 
Since  C02  is  negatively  buoyant,  the  flow  is  highly 
unstable,  and  two-dimensional  visualization  of  the  flow  did 
not  provide  any  evidence  of  the  large-scale  structure  they 
had  predicted. 

Tankin  used  the  same  configuration  as  Ballal  and  Chen 
to  perform  an  extensive  flow  visualization  study.  Air  was 
used  in  the  annulus,  while  N2,  He,  and  C02  were  used 
separately  in  the  jet.  wide  ranges  of  annulus  and  jet  flows 
were  investigated.  TanXin  was  also  able  to  use  a  special 
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technique  that  used  light  scattaring,  Mia  scattering,  to 
determine  tha  shadding  fraquanciaa  of  the  larga  vortex 
structures  evident  in  some  of  tha  flows. 

As  evidenced  by  Tankin' s  work,  velocities  low  enough  to 
avoid  homogeneous  turbulence  produce  large-scale  vortex 
structures  which  persist  far  downstream  of  the  jet's  exit 
plane.  These  flow  patterns  closely  resemble  the  patterns 
seen  in  the  wake  of  a  cylinder  at  low  Reynold**  number  (a  von 
Karman  vortex  street) ,  and  some  of  these  patterns  are  shown 
in  Figure  4.1.  If  the  annulus  velocity  is  fixed,  central 
jet  velocity  increments  result  in  flow  pattern  changes,  some 
of  which  are  quite  remarkable.  These  changes,  called 
bifurcations,  are  the  result  of  interactions  between  the  jet 
and  annulus  flows  and  the  bluff  body. 

To  understand  the  occurrences  in  the  flow  regime  under 
investigation,  flow  visualization,  LDA  measurements,  and 
frequency  information  from  light  scattering  would  be  needed. 
Flow  conditions  at  and  around  the  bifurcation  points  had  to 
be  investigated  quantitatively  as  well  as  qualitatively. 
Extensive  LDA  measurements  could  provide  important 
statistical  quantities  tht  _  would  help  characterize  the 
details  of  the  flowfield.  As  mentioned  earlier,  Ballal  and 
Chen  (1986)  used  the  distribution  of  turbulent  quantities  in 
the  near  field  of  the  slot  jet/bluff  body  flow  to  devise  a 
physical  description  of  the  flow  development.  Knowing  the 
flow  structure  in  the  present  investigation,  one  could  use 
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the  statistical  quantities  to  battar  character! ze  tha  flow 
which  is  depicted  so  wall  in  tha  visualisations. 

Plow  visualisation,  to  include  high  spaed  movies,  would 
ba  dona  not  in  one  plana  but  in  all  throe  planes  of  the  jet. 
Sheet  lighting  the  y-s  (vertical)  plane  (spatial  coordinates 
are  shown  in  Figure  3.8)  shows  the  large-scale  structures  in 
the  flowfield.  Lighting  the  x-i  (lateral)  plane  would 
provide  a  qualitative  estimate  of  the  three-dimensionality 
of  the  flow,  as  would  lighting  the  x-y  (horizontal)  plane. 
Looking  at  the  horisontal  plane  might  also  show  how 
instabilities  in  the  third  dimension  interact  to  cause  or 
assist  in  the  observed  transitions.  Along  these  lines, 
Ghoniam  at  al  (1987)  used  three-dimensional  vortex 
simulation  to  model  the  shear  layer  for  an  axisymmetric  jet. 
They  found  that  the  vortex  rings  formed  in  the  near  field 
became  unstable  (wavy)  to  a  certain  azimuthal  perturbation, 
and  this  caused  the  eventual  breakdown  of  the  vortex  ring. 
The  azimuthal  plane  for  axisymmetric  jets  corresponds  to  the 
horizontal  planw  of  a  planar  jet.  Disturbances  in  the 
horizontal  plane  might  indicate  the  onset  of  an  instability 
that  leads  to  the  bifurcations. 

Mie  scattering  might  provide  information  on  the 
shedding  frequencies  of  the  large-scale  structures  which  may 
provide  important  clues  about  how  this  configuration 
compares  with  other  configurations  in  terms  of  a  parameter 
like  Strouhal  number.  This  comparison  might  also  help 
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explain  the  flow  dynaai.es  and  the  important  factors  in  jat 
development  since  analyses  on  the  simpler  configurations  are 
so  abundant.  If  valid  analogies  can  be  made  for  the  range 
of  jet  velocities  investigated,  one  or  several  pictures  of 
how  the  annulus  and  jet  flows  interact  when  in  the  presence 
of  the  bluff  body  might  emerge.  The  Strouhal  number  for  a 
planar  jet  is  given  by 

f  D 

St  - -  (1.1) 

where  f  is  frequency,  D  is  the  slot  width,  and  is  the  jet 
exit  velocity.  This  is  a  very  important  parameter  for 
categorizing  a  flow  with  frequency  content.  Sato  (1960) 
tested  a  plena  jet  and  observed  that  strouhal  number  was 
independent  of  Reynolds  number  (Re)  for  low  Re  (1500-8000). 
Beavers  and  Wilson  (1970)  used  two-dimensional  and 
axisymmetric  jets  with  Re  from  500  to  3000  and  also  showed 
that  Strouhal  numbers  were  independent  of  Re.  If  Reynolds 
numbers  are  high  enough,  however ,  Strouhal  numbers  for  plane 
jets  (Sato,  1960)  and  axisymmetric  jets  (Becker  and  Massaro, 
1968)  vary  as  the  square  root  of  Reynolds  number.  These 
results  have  been  verified  by  many  other  investigators.  The 
Reynolds  numbers  for  the  present  experiment  (based  on  slot 
width  and  theoretical  (1-D)  jat  exit  velocity)  range  from  35 
to  105.  Although  these  are  quite  low,  one  cannot  be  sure  of 
what  to  expect  since  the  consequences  of  interactions 
between  the  jet  and  annulus  flows  a  re  unknown. 


6 


The  objectives  of  this  investigation  ara  to 
charactarlsa  tha  intaractlons  of  jat  and  annulua  flows  in 
tha  presence  of  a  bluff  body.  With  thaaa  objectives,  a  sat 
of  experiments  was  davisad  to  analyte  tha  flow  pattarn 
changas  and  attaapt  to  dafina  flow  dynamic*  for  this 
two-dimensional  slot  jat/bluff  body  configuration.  Por  tha 
prasant  invast igat ion,  air  was  usad  in  tha  annulus. 
Nitrogan  was  chosen  for  tha  jat  so  buoyancy  would  ba 
nagligibla.  With  a  fixed  annular  valocity,  jat  valocity 
incraaants  that  producad  tha  axpactad  bifurcations  ware 
datarainad  using  a  lasar  sheet  lighting  technique.  Six  jat 
velocities,  with  a  fixad  annulus  valocity,  wara  chosen  for 
investigation.  These  six  conditions  included  jat  velocities 
prior  to,  at,  and  beyond  tha  observed  bifurcation  points. 
Flaw  visualisation  in  two  other  planes  of  the  jet  would  then 
ba  usad  to  assess  tha  threa-diaensi unality  of  tha  jat. 
Lasar  anenoaetry  would  provide  iaportant  statistical 
information  about  tha  flowfields.  Mia  scattering  could  be 
usad  to  identi  j  tha  frequencies  at  which  tha  large 
structures  ara  shad.  This  information  would  ba  coapilad  in 
order  to  provide  a  detailed  picture  of  tha  flowfiald  and  its 
dynaaics.  With  this  plan  of  attack,  than,  tha  coaplate 
character  of  tha  slot  j at/annulus  flow  interaction  aight  be 
batter  understood. 

Evan  with  this  extensive  test  plan  that  includes  some 
sophisticated  techniques,  several  cautions  should  be  noted. 
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Analyses  of  LDA  measurements  of  instantaneous  valocity 
provide  various  statistical  descriptions  of  flow  paraaatars 
(e.g. ,  mean  and  root  aaan  square  velocities,  Reynolds 
stresses,  turbulent  kinetic  energy,  etc) .  Some  of  these 
quantities  will  have  an  entirely  different  interpretation 
for  a  non-turbulant  flow,  and  sose  say  have  no  meaning  at 
all.  Also  important  is  the  fact  that  there  are  two  velocity 
scales  and  two,  possibly  three,  length  scales  for  this 
configuration.  The  jet  and  annulus  velocities  are 
important,  and  ono  sight  imagine  that  the  full  width  of  the 
bluff  body  and  slot,  the  width  of  the  bluff  body  on  only  one 
side  of  the  slot,  and  the  slot  width  could  all  be  important 
lencth  scales.  This  fact  may  make  the  search  for  meaningful 
non-dimensional  parameters,  like  Strouhal  and  Reynolds 
numbers,  more  difficult.  This  certainly  does  not  rule  out 
their  use,  but  these  peculiarities  are  important  and  must 
not  be  ignored. 

In  order  to  boat  describe  the  investigation,  the 
structure  of  this  thesis  will  be  as  fol'.ows.  First,  due  to 
the  unique  nature  of  soma  of  the  experimental  techniques, 
one  chapter  will  be  devoted  to  explaining  the  techniques  and 
their  usefulness.  The  next  chapter  will  bu  used  to  describe 
in  detail  the  experimental  apparatus.  Chaptsr  IV  will 
present  the  background  information  that  motivated  the 
present  investigation;  a  summary  and  discussion  of  previous 
results  obtained  using  the  2-D  centerbody  is  given. 
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Chapter  V  la  a  short  section  describing  the  flow 
conditions  chosen  for  Investigation  and  the  procedures  used 
in  the  exporinents .  Results  of  all  work  will  be  presented 
in  Chapter  VI,  as  will  the  discussion  of  these  results.  The 
conclusions  and  rsoowtndations  will  complete  the  thesis. 


Tna  unique  nature  of  some  of  the  experimental  techniques 
used  in  this  study  demands  that  a  thorough  explanation  be 
given.  Flow  visualisation  which  relies  upon  a  chemical 
reaction  to  produce  subaicron  particles,  a  light  scattering 
technique  to  determine  important  shedding  frequencies,  and 
Laser  Doppler  Anenoaetry  (LDA)  to  determine  mean  and 
fluctuating  velocity  components  were  used  in  the  present 
investigation . 

The  flow  visualisation  technique  has  been  described 
elsewhere  in  detail  (Chen  and  Roquemore,  1986) .  In  this 
section,  the  details  of  the  technique  applicable  to  the 
investigation  will  be  given.  Titanium  tetrachloride 
(TiCl4^)  is  added  to  a  dry  flow,  such  as  the  one  shown  in 
Figure  2.1.  As  the  mixture  exhausts  into  the  ambient  air, 
the  TiCl^  reacts  spontaneously  and  nearly  instantaneously 
with  water  vapor  (H2°(g)^  to  form  T*-°2(»)'  submicron 
particles  which  create  a  light,  white  smoke.  The  chemical 
reaction  is 

TiCl4  +  2H20  4HC1  +  Ti02  (2.1) 

Particles  formed  by  this  reaction  can  then  be  illuminated  by 
various  light  sources  (laser,  strobe,  etc)  to  provide 
precise  visualization.  This  technique  relies  on  having  a 
dry  fluid  in  the  flow  channel  so  preformed  seed  will  not  be 


Figur#  2.1  Flow  Visualization  Tast  Sat-Up 


created.  The  daw  point  temperature  of  tha  fluid  should  be 
checked,  and  tha  valua  should  ba  below  -40°C  (Britton, 
1987) .  It  should  also  ba  noted  that  the  Tio2  particles  are 
near  one  micron  in  size  and  respond  nearly  instantaneously 
to  flow  variations.  The  reaction  given  in  Eg  (2.1)  occurs 
under  nearly  isothermal  conditions,  so  relatively  little 
energy  is  lost  or  gained  by  the  flow.  If  flow  rates  are 
very  low,  however,  this  last  statement  is  invalid.  This 
point  will  be  discussed  in  greater  detail  later. 

Flow  patterns  revealed  by  this  technique  can  be  quite 
interesting.  For  example.  Figure  2.2  shows  some  of  the  flow 
patterns  from  Tankin's  work  (1987).  Illumination  with  a 
sheet  of  laser  light  allows  one  to  photograph  the  flow.  The 
important  feature  of  the  reactive  flow  visualization 
technique  is  that  surfaces  which  have  clear  physical 
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significance  are  illuminated.  Particles  appear  only  at 
those  places  where  the  main  flow  has  contacted  ambient 
fluid.  The  lines  seen  in  the  visualization,  then,  indicate 
the  interfaces  between  jet  and  ambient  fluids. 

The  optical  arrangement  is  shewn  in  Figure  2.3.  A 
Nd:Yag  laser  (Quanta  Ray  Model  DCR-2A)  serves  as  the  light 
source  for  the  visualization.  For  photographs,  such  as  the 
one  shown  in  Figure  2.2,  the  camera  shutter  is  synchronized 
with  a  10  ns  pulse  from  the  laser.  The  pulse  acts  like  a 
flash  for  the  camera.  Photographs  can  also  be  taken  using 
time  exposures  of  flow  fields  " frozen"  by  adjusting  the 
pulse  rate  of  the  laser.  The  sheet  of  light  needed  to 
visualize  the  plane  perpendicular  to  the  slot  jet  is  created 
by  passing  the  laser  output  beam  through  a  telescope  (T) ,  a 
cylindrical  lens  (LI) ,  and  a  magnifying  glass  (MG) .  The 
magnifying  glass  is  used  to  narrow  the  beam  into  a  thin 
sheet.  A  mirror  can  be  used  to  reflect  the  light  sheet  back 
onto  itself,  thereby  increasing  beam  intensity.  To 
visualize  the  plane  parallel  to  that  of  the  jet's  exit 
plane,  one  merely  has  to  rotate  the  cylindrical  lens  90 
degrees . 

Effects  of  the  seed  particles  on  density  and  density 
gradients  had  been  assumed  negligible  (Roquemore,  1987),  but 
this  had  never  been  verified.  A  model,  described  in  detail 
in  Appendix  A,  was  developed  and  verified  the  assumption 
that  the  concentration  of  TiCl4  was  small,  thus  ensuring 
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that  this  fora  of  rsactive  flow  visualization  is 
non-intrusi'/t  (i.a.,  doss  not  altar  the  flow  fiald  that 
would  dsvalop  under  "natural”  circumstances) .  It  should  be 
notad  hara  that  saading  affacts  ara  not  negligible  if  the 
flow  velocities  ara  extremely  small  (  <  5  cm/s) .  Changes  in 
flow  patterns  have  bean  observed  whan  high  seeding 
concentrations  (i.a.,  no  flow  is  being  bypassed  around  the 
TiCl4  seeder)  are  used  in  combination  with  the  lowest 
velocities.  Seed  concentrations  ara  greatest  when  flow 
rates  are  lowest,  and  the  heat  of  reaction  and  density 
gradient  effects  are  the  largest  since  there  is  greater  time 
over  which  the  effects  are  felt.  In  these  instances, 
natural  convection  cannot  be  ignored.  For  this 
investigation,  velocities  greater  than  18.5  cm/s  were 
chosen.  Also,  seed  concentrations  were  minimized  by  setting 
the  valves  on  the  seeder  no  that  as  little  flow  as  possible 
went  to  the  seeder.  These  two  facts  guaranteed  that  no 
pronounced  effects  would  be  seen  for  the  chosen  conditions. 

A  special  type  of  light  scattering,  Mie  scattering,  can 
be  used  to  determine  important  shedding  frequencies.  Large 
scale  structures,  often  known  as  coherent  structures,  are 
seen  in  the  near  field  of  jet  flows  and  are  shed  at  various 
frequencies;  this  method  helps  determine  these  frequencies. 
A  typical  test  arrangement  is  shown  in  Figure  2.4.  The 
reactive  flow  visualization  technique  just  described 
supports  the  scattering  Analysis  by  providing  the  particles 


Figure  2.4  Mie  Scattering  Test  Set-Up 


required  to  scatter  the  laser  light.  The  intersection  of 
the  laser  beam  and  an  imaginary  circle  (one  whose  diameter 
is  that  of  the  fiber  optic  cable)  is  located  4  focal  lengths 
(f)  from  the  end  of  the  optic  cable.  This  intersection 
creates  a  small  spot  which  will  be  called  the  measurement 
spot.  This  spot  must  be  located  somewhere  inside  the 
flowfield  (i.e.,  at  a  location  through  which  Ti02  particles 
will  pass) .  As  the  Ti02  particles  pass  the  measurement 
spot,  the  laser  light  is  scattered.  This  light  is  then 
transmitted  to  a  photomultiplier  via  the  fiber  optic  cable. 
The  signal  from  the  photomultiplier  is  then  amplified  and 
sent  to  either  a  Fast  Fourier  Transform  analyzer  or  a  tape 
recorder,  or  both.  The  tape  recorder  allows  analysis  of  the 
signals  at  a  later  time.  The  resulting  analog  signal  can  be 
analyzed  for  its  frequency  content.  A  similar  method  for 
determining  shedding  frequencies  of  an  axisymmetric  jet  was 
used  by  Becker  and  Massaro  (1968).  The  method,  explained  by 
Becker  et  al  (1967)  used  an  oil  condensation  smoke  to  mark 
the  flow.  Their  scattered-light  technique  used  a  photocell 
and  appropriate  optics  to  generate  a  signal  that  was 
proportional  to  the  instantaneous  point  concentration  of  the 
marking  smoke.  Spectral  analysis  provided  the  shedding 
frequency. 

The  measurement  spot  can  be  located  at  various 
downstream  and  transverse  positions,  so  one  could  determine 
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the  frequencies  at  any  location  in  the  flowfield  where 
particles  pass.  A  sample  time  trace  and  its  corresponding 
frequency  spectrum  are  shown  in  Figure  2.5.  As  evidenced  by 
the  time  trace,  there  is  only  one  meaningful  frequency,  5 
Hr.  The  large  peaks  are  separated  by  approximately  0.2  s, 
as  are  the  smaller  peaks,  and  this  yields  the  frequency  of  5 
Hz.  The  frequency  spectrum  contains  higher  harmonics  only 
as  a  result  of  the  Fourier  transform.  The  time  trace  shown 
in  Figure  2.5(a)  resembles  a  square  ui  j,  and  Fourier 
transforms  of  finite  square  waves  contain  harmonic  frequency 
components  which  are  the  result  of  attempting  to  match  the 
abrupt  changes  at  the  leading  and  trailing  edges  of  the 
square  wave  (Malmstadt  et  al,  1981).  Validation  of  the  Mie 
scattering  method  as  an  effective  means  of  determining 
shedding  frequencies  had  not  been  previously  accomplished. 
Appendix  B  contains  the  details  of  a  set  of  experiments 
designed  to  validate  the  Kie  scattering  technique.  An 
important  factor  to  remember  when  analyzing  Mie  scattering 
information  is  that  interpretation  of  the  results  requires 
the  analog  signal  (time  trace) ,  frequency  spectrum,  and  a 
flow  visualization  photograph.  These  three  pieces  of 
information  are  all  required  to  provide  a  complete  picture 
of  what  is  occurring;  lack  of  one  or  more  pieces  of  this 
information  often  leads  to  unwanted  guess-work. 

Laser  Doppler  Anenometry  (LDA)  was  used  to  determine 
mean  and  fluctuating  components  of  velocity.  The  laser 
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anenometer  is  a  non-intrusive  davica  which  can  provide 
pracisa  quant itativa  measurement  of  valocity.  Measurement 
is  basad  on  the  principle  of  the  Doppler  shift  of  the 
frequency  of  radiation  whan  scattered  from  a  moving  object. 
Small  particles  must  be  distributed  in  the  flow  to  scatter 
the  light.  If  the  particles  are  sufficiently  small ,  the 
slip  velocity  between  the  particle  and  the  fluid  will  be 
small,  and  the  measurement  is  a  valid  measure  of  the  fluid 
velocity.  Optical  arrangements  vary  considerably  based  on 
specific  requirements;  for  example,  see  Holman  (1984).  The 
arrangement  used  in  the  present  investigation  will  be 
described  in  the  following  chapter. 
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The  unique  test  configuration  used  in  the  present 
investigation  is  shown  in  Figures  3.1  and  3.2.  A  slot  (3  mm 
x  150  mm,  aspect  ratio  *50)  is  surrounded  by  a  bluff  body 
(width  of  the  bluff  body  is  6  mm  on  either  side  of  the 
slot).  The  plexiglass  centerbody,  shown  in  greater  detail 
in  Figure  3.2,  is  housed  in  a  large  chimney  section  which 
ends  in  a  24  cm  square  section.  Annulus  flows  are  brought 
around  the  outside  of  the  centerbody  through  the  chimney. 
Two  contraction  sections  are  used  to  minimize  boundary  layer 
effects  in  the  annulus;  contraction  ratio  is  2.68.  A 
square  plexiglass  section  sits  atop  the  chimney  to  minimize 
outside  disturbances. 
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Figure  3.2  Two-Dimensional  Centerbody 
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Call  153  In  Building  490  of  the  Aero  Propulsion 
Laboratory  houses  the  facility  where  all  experiments  ware 
performed.  A  layout  of  the  cell  is  given  in  Figure  3.3. 
Air  was  used  in  the  annulus  and  nitrogen  was  used  in  the 
slot.  The  effect  of  buoyancy  will  be  negligible  with  this 
choice  (various  gases  could  be  used  in  the  central  jet  to 
study  buoyancy  effects,  and  the  facility  can  handle 
combusting  flows;  however,  natural  convection  effects  were 
not  studied) .  Air  could  not  be  used  in  the  central  jet 
since  facility  air  could  not  be  dried  adequately  for  use 
with  the  flow  visualization  and  LDA  seeding.  If  the  fluid 
is  too  moist  and  one  uses  the  TiCl^  for  flow  visualization, 
seed  particles  are  preformed.  This  hinders  clear 
visualization  of  the  flow  structures.  In  addition,  deposits 
may  form  on  the  centerbody  and  annulus  surfaces  and  causes 
the  system  to  be  shut  down  quite  often  for  cleaning.  This 
would  add  an  additional  uncertainty  in  the  results  since  the 
deposit's  effects  are  unknown. 

The  flow  visualization  arrangement  is  shown  in  Figure 
3.4.  The  optical  system  for  visualization  was  described  in 
Chapter  II.  The  annulus  air  is  saturated  with  water  vapor 
to  ensure  the  air  is  humid  enough  for  complete  reaction  of 
the  TiCl4.  Three  Technology  Incorporated  electronic  flow 
control  units  (Model  DMF)  are  used  to  set  and  regulate  the 
required  volumetric  flow  rates  to  within  +  1  percent. 
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The  optical  arrangement  for  the  LDA  io  shown  in  Figure 
3.5.  This  arrangement  has  been  described  in  detail  by 
Ballal  et  al  (1985).  This  LDA  is  able  to  provide 
simultaneous  measurement  of  two  velocity  components.  An  8 
Watt  Argon  ion  laser  (514.5  nm)  is  the  source  for  the  LDA 
system.  Two  measurement  channels  are  separated  by 
polarization.  The  scattered  laser  radiation  is  collected  10 
degrees  off-axis  since  other  systems  on  the  optics  table 
prohibit  using  0  degrees.  The  off-axis  configuration  did 
not  affect  results  (Ballal  et  al,  1986). 

The  data  acquisition  and  processing  system  is  shown  in 
Figure  3.6.  The  counter-processors  analyze  the  signals 


X/2  :  HALF  -  WAVE  PLATE 
8C  :  BRAGG  CELL 
8S  :  BEAM  STEERING  PRISM 

BSh,  BSv  •  BEAM  SPLITTER  (HORIZONTAL,  VERTICAL) 
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(o)  SCHEMATIC  DIAGRAM  OF  THE  LDA  OPTICS  ARRANGEMENT 
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(b)  SCHEMATIC  OF  RHOMB  PRISM  SET  SHOWING  INCOMING  BEAMS 
(POSITION  I),  OUTGOING  BEAMS  (POSITION  2)  AND  THEIR 
RESPECTIVE  CENTERS. 


Figure  3.5  Two -Component  LDA  System  (Optics) 
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Figure  3.6  Data  Acquisition  and  Processing  System 


provided  by  each  photodetector  (one  for  each  of  the  two 
channels)  of  the  LDA.  The  Doppler  signal  from  the 
photodetector  is  amplified,  filtered,  and  then  converted 
into  a  pulse  for  digital  processing  by  the 
counter-processors.  These  units  process  only  those  signals 
that  have  exceeded  some  threshold,  thus  ensuring  good  signal 
quality.  The  processor  requires  32  cycles  in  an  LDA  burst 
(a  burst  occurs  when  a  particle  passes  through  the 
measurement  volume).  The  Bragg  cell  shown  in  Figure  3.5 
helps  meet  this  criterion.  The  cell  operates  at  40  MHz, 
thus  shifting  the  zero  velocity  signal  to  40  MHz.  The 
Doppler  signal  lies  on  one  side  or  the  other  of  this 
frequency,  depending  on  the  scatterer's  direction.  The 
frequency  shift  also  increases  the  number  of  cycles 
available  for  analysis.  For  example,  5  m/s  velocities  would 
yield  Doppler  bursts  with  a  maximum  of  10  /is.  A  typical 
burst  would  be  3  /is;  at  40  MHz,  120  cycles  would  be  provided 
for  analysis.  This  set-up  provides  poor  resolution  for  low 
speed  flows  (<1  m/s) ,  so  a  phase-locked  down-mixer 
electronically  shifts  the  signal  from  40  MHz  to  5  MHz  or  10 
MHz  (user  selectable).  At  10  MHz,  velocity  resolutions  of 
1.1  cm/s  vertically  and  1.4  cm/s  horizontally  are  obtained. 

The  processed  signals  are  then  sent  to  the  LDA 
interfaces  which  record  arrival  time  and  the  number  of 
cycles  in  the  burst.  The  counter-processor  and  clock  data 
are  sent  to  a  mini-computer  where,  when  4096  data  points  at 
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a  specific  location  have  been  recorded,  a  preliminary 
analysis  is  done  which  provides  estimates  of  mean  velocity 
and  the  root  mean  square  (rms)  of  the  fluctuating 
velocities.  Data  are  then  saved  on  magnetic  tape,  and  the 
system  moves  to  the  next  point  and  begins  to  take  data 
again.  This  process  is  all  computer-controlled. 

The  chimney  and  centerbody  are  located  on  a  traversing 
table  which  is  shown  in  Figure  3.7.  Vertical  traverse  of  2 
feet  and  horizontal  traverses  of  9  inches  (both  ways)  are 
possible.  By  moving  the  burner  and  not  the  optics, 
complexity  of  the  optics  is  minimized  and  repeated 
realignments  and  inaccuracies  are  avoided.  The  traverse 
mechanism  controls  position  within  ±  10  jim.  The  spatial 
coordinate  system  is  shown  in  Figure  3.8,  and  the 
corresponding  coordinates  for  velocity  are  given  in  Figure 
3.9.  These  systems  should  be  noted  well.  Also,  the  terms 
transverse  and  radial  will  be  used  interchangeably.  This 
notational  inconvenience  may  arise  at  a  later  time,  but  all 
attempts  will  be  made  to  ensure  clarity  in  data  presentation 
ana  discussions. 

The  seeding  used  for  the  LDA  deserves  special 
attention.  Seeders  must  produce  a  steady  flux  of  micron  or 
submicron  size  particles  and  must  be  controllable.  Adequate 
and  steady  seeding  of  proper  size  particles  ensures  good 
signal  quality  and  accurate  velocity  measurement.  High 
seeding  rates  must  be  avoided  since  multiple  particles  will 
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Figure  3.7  Three-Axis  Burner  Traversing  Platform 
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Figure  3.8  Spatial  Coordinate  System 
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scatter  radiation  in  diffarant  fashions  and  yiald  unusable 
rasults.  Low  saading  ratas  rasult  in  data  acquisition  which 
is  vary  slow,  so  tins  to  acquira  the  desired  information  can 
become  vary  large.  Spectral  information  is  also  often 
important  in  analyses,  and  low  saading  gives  data  rates 
inadequate  to  gat  frequency  or  power  spectra.  These 
problems  are  amplified  whan  low  velocity  flows  are  under 
study. 

The  test  arrangement  used  for  LDA  measurements  is  shown 
in  Figure  3.10.  The  high  output  fluidized  bed  seeder  was 
used  for  the  annulus  flow  and  provided  consistent  seeding. 
The  major  drawback  of  many  fluidized  bed  seeders  is  their 
inability  to  obtain  the  desired  particle  generation  rate. 


Figure  3.10  LDA  Flow  System  Arrangement 
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This  experiment  was  dons  at  a  constant  annular  velocity,  so 
after  the  initial  setting,  there  were  no  more  adjustments 
and  the  problem  was  avoided.  Tha  low  output  seeder  used  for 
the  central  jet  was  a  significant  problem  for  the  LDA 
measurements.  Very  low  volumetric  rates  (-  5-15  slpm)  made 
it  nearly  impossible  to  get  adequate  seeding  rates.  At 
first,  the  pipe  shown  in  Figure  3.11  was  used.  Some  of  the 
jet  flow  could  be  bypassed  around  the  seeder  if  necessary. 
Frequent  clogging  of  the  screens  soon  prohibited  its  use.  A 
cyclone-type  seeder  was  then  created  by  Mr  Ron  Britton  of 
the  Aero  Propulsion  Laboratory.  This  device,  shown  in 
Figure  3.12,  uses  a  motor  to  stir  a  bed  of  A1203  particles. 
Air  is  brought  in  the  seeder  so  a  cyclone  is  created  which 
will  take  with  it  the  aluminum  oxide  particles.  Flow  could 
also  be  bypassed  around  the  seeder  to  whatever  degree 
desired.  This  configuration  worked  extremely  well  and 
provided  excellent  seeding  over  the  range  of  flows 
investigated . 

When  data  have  been  collected,  various  computer 
programs  are  available  to  analyze  the  data.  These  programs 
have  been  developed  by  personnel  from  the  University  of 
Dayton  Research  Institute  (UDRI) .  The  user  has  great 
flexibility  since  all  the  programs  are  menu  driven.  Various 
averaging  schemes  and/or  data  smoothing  can  be  used,  and 
statistical  descriptions  of  the  data  such  as  mean  velocity, 
root  mean  square  (rms)  of  the  fluctuating  velocity, 
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Figure  3.12  Low  Output  Cyclone  Seeder 

covariance ,  turbulent  kinetic  energy ,  skewness,  kurtosis, 
probability  density  functions,  and  spectral  information  ctin 
be  obtained.  Techniques  for  data  analysis  and  presentation 
wi  be  described  as  appropriate  throughout  the  remainder  of 
this  text. 
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IY*.  Preliminary  Investigations 


Motivation  for  the  present  investigation  cane  from  the 
work  of  Tankin  (1987).  The  centerbody  previously  described 
was  used  in  Tankin' s  work,  and  the  primary  emphasis  was  on 
flow  visualization.  Only  a  small  portion  of  quantitative 
data  was  taken.  Various  annulus  velocities  were  run,  and 
jet  gases  of  He,  N2,  and  C02  were  run  over  a  broad  range  of 
velocities.  A  typical  sequence  of  photographs  is  shown  in 
Figure  4.1.  These  photographs  show  N2  in  the  central  jet. 
After  the  annulus  velocity  is  set  at  its  prescribed  value, 
the  jet  flow  is  increased  from  its  minimum  value  until  all 
jet  flows  have  been  run.  Jet  to  annulus  velocity  ratios 
typically  varied  from  0.25  up  to  5.0  or  greater  to  ensure 
wide  coverage  of  flow  conditions.  Photographs  were  taken  at 
each  condition  using  the  laser  sheet  lighting  technique 
described  in  Chapter  II.  Theoretical  (1-D)  exit  velocities 
for  the  jet  and  annulus,  those  calculated  using 

Q 

U  -  -  (4.1) 

A 

where  Q  is  volumetric  flow  rate  and  A  is  area,  will  be 

denoted  by  and  U  ,  respectively, 
j  a 

Several  key  features  of  the  flow  visualization 
technique  should  be  pointed  out  here.  As  discussed  in 
Chapter  II,  the  areas  where  the  Ti02  particles  appear  mark 
the  boundaries  between  jet  and  annulus  fluid.  Several 


Figure  4.1  Flow  Visualization  from  Tankin  (1987) 


Figure  4.1  Flow 


important  terms  will  be  used  when  analyzing  the  flow 
visualization  photographs,  and  these  are  provided  in  Figure 
4.2.  One  can  clearly  see  that  vortex  orientation  is  "down" 
for  Uj  >  37.0  cm/s,  while  it  is  "up"  for  <  37.0  cm/s. 
The  "trunk"  varies  considerably  for  the  eleven  cases.  When 
Uj  «  37.0  cm/s,  a  trunk  does  not  even  exist,  only  interface 
lines.  Familiarity  with  these  terms  will  be  important  in 
later  analyses  of  photographs. 

The  large  scale  structures  in  Figure  4 . 1  are  quite 
interesting.  More  intriguing,  however,  are  the  drastic  flow 
pattern  changes  which  are  seen.  The  cases  where  Uj  =  22.2 
cm/s  and  »  37.0  cm/s  are  particularly  interesting.  At 
■  22.0  cm/s,  the  trunk  has  become  quite  large,  and  there 

appears  to  be  a  stagnation  point  located  at  the  base  of  the 
trunk  (this  apparent  stagnation  was  identified  when  viewing 
the  high  speed  movies).  At  Uj  ■  37.0  cm/s,  almost  parallel 
interface  lines  are  produced.  Small  perturbations  (ones 
that  are  obviously  not  amplified)  prevent  the  lines  from 
being  perfectly  straight. 

The  two  flow  patterns  just  discussed  are  very  sensitive 
to  changes  in  jet  velocity.  The  patterns  can  change 
drastically  with  only  a  small  change  in  jet  velocity.  Once 
the  pattern  is  established,  however,  it  is  very  stable.  The 
pattern  can  be  perfectly  frozen  using  the  variable  pulse 
rate  of  the  laser  or  a  variable  frequency  strobe.  The  laser 
provides  a  thin,  two-dimensional  sheet  of  light,  while  the 


38 


lfc>  Down 

Vortex  Orientation 


Figure  4.2  Flow  Visualization  Terminology 
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strobe  can  be  used  to  illuminate  the  flow  from  a  broad 
perspective.  These  strong  points  can  be  used  in  combination 
to  provide  a  detailed  picture  of  the  flow  in  all  dimensions. 
The  pulse  rate  of  the  laser  can  be  selected  in  a  small  ranqe 
of  frequencies  (these  vary  depending  on  the  desired  pulse 
strength  which  is  in  Joules  per  pulse)  which  are  very  close 
to  the  frequencies  at  which  the  observed  structures  are 
shed;  frequency  of  the  strobe  could  be  varied  from  1.67  to 
416.67  Hz.  Varying  the  source's  pulse  rate  results  in 
structures  which  appear  to  be  going  forward  or  backward  or 
can  be  exactly  frozen.  The  apparent  velocity  of  the 
structures  depends  on  the  difference  between  the  source  and 
shedding  frequencies.  If  shedding  is  very  repetitive,  the 
flow  can  be  perfectly  frozen  by  proper  selection  of  the 
light  source  frequency.  The  ability  to  freeze  a  certain 
pattern,  then,  suggests  a  particular  stability  of  the  flow. 

The  other  flow  patterns  of  Figure  4 . 1  change  as  flow 
velocities  change,  but  the  basic  patterns  remain  the  same. 
For  example,  the  cases  where  Uj  *  27.8  and  33.3  cm/s  show 
patterns  which  are  not  identical,  but  certainly 
qualitatively  similar.  Using  the  terminology  of  Perry  and 
Chong  (1987),  patterns  are  considered  topologically 
equivalent  if  they  can  be  distorted  into  one  another  by  a 
stretching  process.  If  we  picture  each  pattern  on  a  rubber 
sheet,  then  the  first  four  patterns  could  be  created  by 
arbitrarily  stretching  the  first  pattern.  The  fifth 
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pattern,  however,  could  not  be  created  by  stretching  of  any 
of  the  other  patterns.  This  drastic  flow  pattern  change,  a 
bifurcation,  is  indicative  of  a  structurally  unstable 
critical  point  (Perry  and  Chong,  1987) .  Another  bifurcation 
seems  to  have  occurred  when  Uj  =  37  cm/s.  It  was  also  noted 
by  Perry  and  Chong  (1987)  that  near  a  structurally  unstable 
point,  a  slight  change  in  one  of  the  relevant  parameters 
could  cause  a  spectacular  change  in  the  way  a  jet  or  wake 
distributes  and  transports  (transfers)  its  mass,  momentum, 
and  energy.  Although  the  work  of  Perry  et  al  has  not 
focused  on  configurations  like  the  one  under  investigation, 
it  was  felt  that  the  physics  of  the  flows  might  be  similar 
enough  to  provide  at  least  some  insight  into  flowfield 
dynamics.  Interpretation  in  terms  of  bifurcations  may  help 
in  understanding  the  interaction  of  the  bluff  body  and  jet 
flows.  Unfortunately,  no  quantitative  analyses  could  be 
done  since  they  are  extremely  complex. 

Having  established  the  motivation  for  the  present 
experiment,  flow  conditions  and  experimental  procedure  had 
to  be  specified.  Since  flow  visualization,  LDA,  and  Hie 
scattering  would  be  used,  conditions  had  to  be  chosen 
judiciously  because  of  time  constraints.  These  factors  and 
the  experimental  procedures  are  discussed  in  the  next 
chapter. 
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3L.  Test  Conditions  and  Experimental  Procedure 


To  accomplish  the  objectives  stated  in  Chapter  I,  a 
test  plan  was  selected  that  would  provide  adequate 
characterization  of  the  flowfield  interactions.  The 
limiting  factor  in  choosing  the  number  of  cases  to 
investigate  was  the  LDA  measurements.  Because  of  low 
seeding  rates,  a  test  plan  had  to  be  created  that  would  be 
coarse  enough  to  ensure  deadlines  were  met  but  fine  enough 
to  guarantee  adequate  coverage  of  the  flow.  Six  conditions 
were  chosen  for  investigation  and  are  given  in  Table  I. 
These  captured  both  bifurcation  points  as  well  as  conditions 
before  and  after  these  points. 


TABLE  I  Flow  Conditions 


Condition 

Ua(cm/s) 

U . (cm/s) 

(Qa/Qj)  slpm 

1 

27.5 

18.5 

(330/5) 

2 

27.5 

22.2 

(330/6) 

3 

27.5 

27.8 

(330/7.5) 

4 

27.5 

33.3 

(330/9) 

5 

27.5 

37.0 

(330/10) 

6 

27.5 

51.8 

(330/14) 

It  was  decided  to  repeat  six  of  Professor  Tankin' s 
conditions,  partly  because  of  the  familiarity  with  his  work. 
Different  annular  velocities  were  tried  in  preliminary 
visualization  studies,  but  jet  flow  variations  led  to  the 
same  sequence  of  patterns  observed  earlier.  A  given 
pattern,  however,  would  not  repeat  itself  at  a  specified 
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velocity 

behavior 


This 


ratio  U./U  for  different  values  of  U  . 

J  a  a 

was  caused  by  the  turbulence  intensity  in  the 
annulus.  The  turbulence  intensity  for  the  annulus  was 
directly  proportional  to  the  volumetric  flow  rate. 
Therefore,  even  though  the  patterns  were  highly  repeatable, 
they  were  also  very  susceptible  to  variations  in  annulus 
turbulence  levels.  For  the  2-D  centerbody  configuration, 
complete  break-up  of  structures  was  seen  at  approximately  Q. 

a 

*  500  slpm  (the  annulus  had  become  "fully  turbulent"  at  this 
point) .  Although  one  would  like  as  high  a  velocity  as 
possible  to  ensure  good  data  rates  for  the  LDA,  it  was 
decided  to  avoid  high  volumetric  flow  rates  in  the  annulus. 
Professor  Tankin' s  conditions  offered  good  middle  ground. 

Another  factor  may  contribute  to  the  discrepancy  that  a 
specific  velocity  ratio  does  not  produce  the  same  flow 
pattern  when  Ua  is  changed.  For  a  fixed  Uj/Ua,  must 
change  when  U&  is  changed.  At  low  velocities,  one  might 
expect  a  parabolic  jet  exit  profile,  while  high  jet 
velocities  should  tend  toward  a  top-hat  profile.  Therefore, 
it  may  be  important  to  look  at  exit  mean  velocity  profiles 
in  order  to  ascertain  the  correct  exit  velocity. 

The  repeatability  of  the  sequence  of  flow  patterns  has 
been  mentioned  several  times  and  deserves  a  short  comment 
here.  The  qualitative  evolution  of  these  patterns  at 
different  annular  velocities  is  nearly  identical. 
Quantitative  results  will  differ,  but  the  test  conditions 
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are  representative  of  a  general  evolution  of  the  flow. 

Studying  one  set  of  conditions  instead  of  another,  then, 

should  not  lead  to  different  characterizations  of  the 

flowfield  dynamics  —  an  encouraging  result.  It  should  also 

be  noted  here  that  the  repeatability  of  the  flow  patterns 

for  a  particular  TJ.  (with  a  fixed  U  )  is  perfect.  In  other 

J  a 

words,  for  a  fixed  annular  velocity,  if  one  obtains  a 
certain  pattern  at  a  Uj,  changes  flow  conditions,  then 
returns  to  the  original  ,  the  original  pattern  will  always 
return . 

Once  the  flow  conditions  were  determined,  the  LDA  data 
acquisition  grid  was  chosen.  It  may  be  appropriate  at  this 
time  to  review  Figure  3 . 8  since  it  contains  the  spatial 
coordinate  system.  Only  one  lateral  position,  the  center  of 
the  slot,  was  used.  At  a  given  axial  position,  a  "run" 
consists  of  a  transverse  (Y  direction)  scan.  The  Y 
positions  chosen  for  each  run  are  summarized  in  Table  II. 


TABLE  II  Transverse  Scan  Locations 


Y  RANGE  (mm) 

INCREMENT  4?  (mm) 

(-30)  to  (-15) 

5 

(-15)  to  (-9) 

2 

(-9)  to  (+9) 

1 

(+9)  to  (+15) 

2 

The  edge  of  the  annulus  extends  to  -42  min,  but  it  was  shown 
by  Tankin  that  edge  effects  for  |Y|  >  30  mm  did  not  impact 
the  flow  near  the  bluff  body.  The  contraction  section  was 
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effective  in  providing  a  top-hat  velocity  profile  at  the 
annulus  edge.  Due  to  the  symmetry  of  the  configuration,  a 
full  transverse  scan  was  not  necessary.  Each  scan  consisted 
of  26  points  where  the  u  and  v  (axial  and  transverse, 
respectively)  velocities  would  be  measured.  Lateral 
velocity  could  not  be  measured  since  the  LDA  system  could 
only  measure  two  velocity  components.  Axial  positions  where 
the  transverse  scans  would  be  taken  were  chosen  at  z  **  6, 
10,  20,  40,  and  60  mm.  A  true  exit  profile  (z«0)  could  not 
be  obtained  because  the  centerbody  interfered  with  the  LDA 
beams  below  6  mm.  The  combination  of  transverse  scans  at  5 
axial  positions  resulted  in  measurements  at  130  points  in 
the  flowfield.  This  provided  coverage  dense  enough  to 
indicate  how  quantities  are  changing.  Because  of  the  lack 
of  structures  at  low  z  for  Uj  ■  51.8  cm/s  (see  Figure  4.1), 
data  were  taken  at  z  -  10,  20,  60,  and  120  mm  for  this  case 
only. 

Flow  visualization  in  all  three  planes  of  the  jet  was 
accomplished  for  the  conditions  given  in  Table  I.  For  the 
vertical  sheet  lighting,  photographs  at  each  of  the  six 
conditions  were  taken.  This  was  also  done  for  the  lateral 
sheet  lighting.  For  the  horizontal  sheet  lighting, 
photographs  were  taken  at  each  of  the  z  locations  specified 
earlier,  thus  providing  a  picture  of  the  influence  of  the 
third  dimension  as  axial  position  increased.  The 
photographic  summary  is  provided  in  Table  III. 
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TABLE  III  Photographic  Summary 
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The  Mi®  Scattering  Method  (MSM)  was  used  to  determine 
the  shedding  frequencies  of  the  large-scale  structures.  The 
optics  for  the  LDA  were  used  instead  of  the  arrangement 
described  in  Chapter  II  (which  was  only  used  to  validate 
this  technique) .  The  technique  itself  is  unchanged. 
Stroboscopic  lighting  was  first  used  to  ensure  the 
measurement  spot  was  in  a  region  where  the  vortices  would 
pass.  Once  in  this  region,  the  centerbody  traverse  was 
moved  until  maximum  signal  strength  was  observed  on  the 
oscilloscope  (which  was  part  of  the  LDA  arrangement) .  The 
high  voltage  source  was  then  adjusted  to  provide  a  signal 
having  a  1.5  volt  peak  which  would  be  sent  to  the  FM  tape 
recorder.  Data  were  recorded  at  axial  positions  of  20,  40, 
and  60  mm  for  the  first  five  jet  velocities  (up  to  37.0 
cm/s) ,  and  at  40,  60,  and  120  mm  for  the  case  where  « 
51.8  cm/s.  This  distribution  of  stations  was  chosen  to 
coincide  with  LDA  measurement  locations,  but  measurements  at 
lower  z  locations  (z  <  20  mm)  were  deleted  since  information 
about  the  large  structures  is  desired,  and  these  are  not 
evident  very  near  the  jet. 

A  case  that  illustrates  the  possible  interpretation 
errors  in  using  frequency  analysis  will  be  discussed  here. 
Figure  5.1  shows  the  time  trace  and  its  associated  frequency 
spectrum  at  z  *  40  mm  for  the  case  where  -  18.5  cm/s. 
The  shedding  frequency,  while  it  appears  from  the  frequency 
spectrum  to  be  12.25  Hz,  is  actually  6.125  Hz.  Looking  at 
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Volts 


(b)  Frequency  Spectrum  of  Signal  in  (a) 

Figure  5.1  Potential  Misinterpretation  of  Mie 
Scattering  Results 


the  tin*  trace  in  Figura  5.1(a),  one  can  see  that  a  full 
waveform  (which  in  this  example  consists  of  two  peaks  of  two 
different  amplitudes)  repeats  itself  every  0.16  seconds, 
thus  a  frequency  of  6.125  Hz.  The  Fourier  analyzer, 
however,  has  seen  the  flat  zones  at  the  top  of  the  signal 
(two  per  period)  and  correlated  them  so  that  they  have  a 
frequency  of  12.25  Hz.  This  correlation  is  a  strong  one, 
and  the  frequency  spectrum  indicates  this.  For  a  full 
discussion  of  Fourier  transforms  of  finite  signals,  see 
Malmstadt  et  al  (1981).  In  most  other  cases,  the  Mie 
scattering  data  required  proper  interpretation  to  determine 
the  shedding  frequency.  This  emphasizes  the  need  to  examine 
the  signal's  time  trace.  Photographs  were  also  important  in 
resolving  this  apparent  discrepancy.  Measuring  the  spacing 
of  the  structures  passing  at  40  mm,  a  shedding  frequency  of 
12.25  Hz  would  imply  that  the  convecting  velocity  of  the 
structure  was  greater  than  both  the  annulus  and  jet 
velocities.  Using  6.125  Hz  gives  a  result  consistent  with 
the  photograph  and  the  mean  axial  velocity. 
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YLl  Results  and  Discussion 


Flow  Visualization 

Flow  visualisation  in  the  three  planes  of  the  jet  were 
done.  Since  the  photographs  taken  in  the  y-z  plane  (those 
that  show  the  large-scale  structures)  were  expected  to  be 
very  similar  to  those  of  Tankin,  it  is  more  interesting  to 
view  visualizations  in  the  other  two  planes  first.  The 
lateral  (x-z)  plane  visualizations  are  shown  in  Figure  6.1. 
These  photographs  identify  a  problem  with  the  centerbody 
which  had  been  noted  previously  under  stroboscopic  lighting. 
The  flow  is  not  two-dimensional  (i.e.,  there  is  not  a 
top-hat  profile  in  the  lateral  direction) ,  particularly  at 
low  velocities.  When  the  problem  was  first  noted,  it  was 
felt  that  the  flow  dividers  (see  Figure  3.2)  were  not 
equispaced,  thus  distributing  the  flow  improperly.  This 
problem  was  fixed  by  realigning  the  dividers  at  the  bottom 
of  the  centerbody.  When  reinstalled,  the  flow  patterns  of 
Figure  6.1  were  unchanged.  The  problem  may  be  one  of  flow 
diffusion.  The  area  ratio  of  each  of  the  4  flow  segments  is 
greater  than  4.0.  At  very  low  velocities,  low  energies,  the 
flow  cannot  diffuse  efficiently.  The  honeycomb  and  screens 
in  the  chimney  section  should  remove  any  lateral  swirl 
present  in  the  annulus  flow,  but  this  too  could  be 
contributing  to  the  problem.  This  problem  was  not 
investigated  further.  It  is  interesting  to  note  that  the 
flow  for  the  worst  case,  =  18.5  cm/s,  still  maintained  a 
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Figure  6.1  Lateral  Plane  Flow  Visualization 


Figure  6.1  Lateral  Plane  Row  Visualization 


Figure  6.1  Lateral  Plane  Flow  Visualization 


two-dimensionality  (away  from  the  edges) ,  just  that  it 
occurred  in  a  plane  that  now  had  the  slope  of  the  lines 
shown  in  Figure  6.1.  Using  the  strobe  as  a  light  source, 
one  could  look  down  the  centerline  of  the  vortex  sheet 
(viewing  from  the  left  looking  at  Figure  6.1)  at  the 
appropriate  angle,  and  it  was  evident  the  flow  was 
two-dimensional.  It  is  best  to  use  a  term  like  "quasi-2-D" 
or  "piecewise  2-D"  to  describe  the  flows  in  this  experiment. 
In  their  study  of  a  free  shear  layer,  Lasheras  et  al  (1986) 
found  that  the  characteristics  of  the  Kelvin-Helmholtz 
roll-ups  in  the  shear  layer  were  virtually  unaffected  by  the 
presence  of  three-dimensionality.  Although  true 
two-dimensionality  may  never  be  achieved  for  this 
configuration,  all  attempts  were  made  to  minimize  spurious 
contributions  to  the  three-dimensional  behavior.  It  was 
felt  that  this  had  been  at  least  roughly  achieved,  so  the 
experiments  proceeded  using  the  existing  configuration  and 
test  conditions. 

Because  of  the  three-dimensionality  (tilt)  of  the  flow, 
the  photographs  in  the  horizontal  plane  were  not  able  to 
provide  much  insight  to  how  instabilities  in  that  plane 
affected  the  transitions.  The  photographs  at  the  various  z 
locations  for  Uj  =18.5  cm/s  and  37.0  cm/s  are  given  in 
Figures  6.2  and  6.3,  respectively.  For  =  18.5  cm/s,  it 
is  nearly  impossible  to  tell  if  the  observed  pattern  is  a 
result  of  a  true  instability  and  vortex  roll-up  (analogous 
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(e)  z-60  mm 


■ 

(d)  z  «  40  mm 

(a)  z  -  6  mm 

Figure  6.2  Horizontal  Flow  Visualization 
(Uj  *  18.5  cm/s) 
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to  a  Tollmein-Schlichting  wave) ,  or  the  result  of  fluid 
motion  in  the  lateral  direction  induced  by  the  experimental 
imperfections.  Unfortunately,  it  seems  to  be  the  latter. 
Viewing  Figure  6.3  does  show,  as  did  Figure  6.1,  that  the 
jet  is  well-behaved  and  two-dimensional  at  the  high  flow 
rates  (even  as  high  as  z  =60  mm)  .  Due  to  the 
three-dimensional  nature  of  the  flow,  however,  specific 
analyses  would  be  pure  speculation  at  this  time  and  are 
avoided. 

Some  general  comments  on  these  visualizations  are 
appropriate.  As  alluded  to  earlier,  one  can  see  how  far 
downstream  the  two-dimensionality  of  the  flow  prevails.  2-D 
flow  would  imply  that  we  would  see  a  rectangle  in  the 
horizontal  plane.  As  one  can  see  in  Figure  6.2,  flow  at  Uj 
=  18.5  cm/s  already  has  developed  a  small  amount  of 
waviness  at  z  =  6  mm.  The  same  is  true,  but  to  a  lesser 
degree,  for  jet  velocities  of  22.2  and  27.8  cm/s  (these 
photographs  are  not  shown  here) .  Higher  velocities  result 
in  "two-dimensional”  flows  at  this  z  location.  As  one  goes 
to  greater  and  greater  axial  positions,  the  influence  of  the 
third  dimension  becomes  larger  and  larger. 

An  interesting  behavior  is  observed  in  all  photographs 
of  the  horizontal  plane.  At  the  ends  of  the  slot,  a 
rounding  has  occurred  which  makes  the  shape  of  the 
illuminated  particles  much  like  a  dog's  bone.  One  might 
think  about  this  effect  in  the  following  way.  Imagine 
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starting  with  an  axisymmetric  jet.  The  pattern  seen  in  its 
horizontal  plane  would  be  a  circle.  Distortion  of  the  jet 
by  squeezing  its  two  sides  would  begin  to  produce  ovals 
which  would  become  more  and  more  oblong  as  the 
squeezing  continued.  At  some  point  a  rectangle  (similar  to 
the  plane  jet)  could  be  achieved.  At  the  ends,  some 
rounding  would  still  exist,  but  these  effects  will  have 
been  confined  to  a  small  region  if  squeezing  was  stopped  at 
a  high  aspect  ratio.  The  analogy  here  is  only  meant  as  a 
tool;  this  picture  of  the  plane  jet  may  be  useful  for 
conceptualizing  various  aspects  of  the  flow. 

The  vertical  sheet  photographs  are  given  in  Figure  6.4. 
These  compare  quite  well  with  Tankin' s  work  (see  Figure 
4.1),  as  they  should.  The  case  where  uj  “  37.0  cm/s  is 
where  a  bifurcation  occurs,  and  the  new  test  conditions 
produced  the  pattern  seen  in  Tankin' s  photograph  at  a 
slightly  different  jet  velocity  (<2  cm/s  difference).  There 
were  subtle  differences  between  the  present  configuration 
and  the  one  used  by  Tankin,  and  this  would  explain  the 
observed  differences  in  flow  patterns.  The  contraction 
sections  used  in  the  two  experiments  were  different.  The 
contraction  ratio  for  Tankin' s  work  was  4.5,  and  it  was  2.68 
in  the  present  work.  This  difference  meant  the  annular  flow 
rate  Qa  had  to  be  increased  to  meet  the  requirement  that 
velocities  in  the  two  experiments  were  the  same.  As 
discussed  earlier,  this  will  change  the  turbulence  intensity 
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Plan*  Flow  Visualization 


in  the  annulus.  Even  with  these  differences,  the  degree  to 
which  the  flows  are  similar  indicates  that  a  similar 
qualitative  evolution  of  these  flowfields  is  occurring. 

llBft  MAaauraniftnta 

The  LDA  system  described  earlier  was  used  to  gather 
instantaneous  axial  and  transverse  velocity  data.  With  the 
data  obtained,  one  could  evaluate  mean  velocities,  the  root 
mean  square  of  the  velocity  fluctuations,  normalized 
Reynolds  stress  (covariance) ,  turbulent  kinetic  energy, 
kurtosis,  and  skewness.  One  could  also  get  the  probability 
density  function  (PDF)  and  the  spectral  information  from 
data  obtained.  Many  methods  of  data  presentation  were 
available. 

The  first  problem  which  is  evident  is  that,  while  the 
flow  is  clearly  non-turbulent,  only  quantities  typically 
used  to  analyze  turbulent  flows  are  available.  Except  for 
mean  velocities  and  spectral  analyses,  these  quantities 
demand  new  interpretation  if  they  will  be  used  to  help 
determine  the  character  of  the  flow. 

Mean  values  have  a  straightforward  interpretation.  The 
present  significance  of  root  mean  square  velocities  can  be 
analyzed  by  considering  the  following  idealized  case. 
Consider  the  transverse  velocity  for  a  2-D  vortex  (assumed 
to  be  a  solid  body)  translating  in  only  the  axial  direction. 
If  the  measurement  volume  splits  the  vortex  as  it 
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translates,  ona  might  expect  the  instantaneous  velocity 

signal  to  look  like  a  sine  wave.  The  concept  is  shewn  for 

two  vortices  in  Figure  6.5.  The  tilda  denotes  an 

instantaneous  value.  At  tQ,  the  first  vortex  passes  the 

measurement  spot  with  its  minimum  velocity.  Velocity  then 

increases  through  zero  (center  of  vortex)  until  the  bottom 

of  the  first  vortex  is  reached.  The  space  between  vortices 

is  assumed  to  produce  a  smooth  velocity  gradient  when  going 

from  the  maximum  (end  of  the  first  vortex;  back  to  the 

minimum  (beginning  of  the  second  vortex)  velocity.  The 

process  then  repeats  itself  for  the  second  vortex.  Since 

v _  is  obtained  from  the  data  analysis,  and  a  functional 

rms 

form  for  has  been  assumed,  v.„  „  can  be  related  to  the 

rms 

rotation  speeds  of  the  vortices.  For  a  sine  wave,  the 
maximum  amplitude  and  its  root  mean  square  are  related. 
Using  the  analogy  to  the  rigid  rotator,  one  can  say 


max' 


VT 


rms 


(6.1) 


Flow  visualization  can  provide  vortex  size,  r.  Angular 
speed,  u>,  could  be  estimated  using 


< t) 


|v. 


max' 


r 


(6.2) 


A  non-zero  mean  velocity  will  affect  this  calculation 
slightly  since  this  means  the  solid  body  is  also  translating 
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Figure  6.5  idealized  Instantaneous  Transverse 
Velocity  for  Two  Solid  Body 
Rotators  Translating  Uniformly 


in  the  Y  direction.  For  this  analysis,  however,  the  effect 
will  be  ignored.  This  method  could  provide  insight  into  the 
dynamics  of  the  vortices. 

The  Reynold  stress,  u'v' ,  normalized  by  the  product 
2  o 

(u*  •  v'  ),  the  variances  of  u  and  v,  is  called  the 
covariance.  It  is  a  value  between  0  and  +1  that  identifies 

the  degree  to  which  u _  and  v _  are  correlated.  If  a  and 

b  represent  two  signals,  if  a  is  always  positive  when  b  is 
negative,  then  the  parameters  are  perfectly  negatively 
correlated  (i.e.  covariance  is  -1).  If  the  signals  are  both 
positive  most  of  the  time,  the  signals  are  said  to  be 
strongly  positively  correlated  (i.e.  covariance  greater 
than  zero  and  near  unity) .  Two  completely  random  signals 
should  yield  a  covariance  of  zero.  For  non-turbulent  flows, 
it  will  be  important  to  look  at  covariance  as  a  correlation 
parameter  only,  not  as  a  Reynolds  stress. 

with  these  ideas  in  mind,  the  data  are  presented  in 
Figures  6.6  through  6.10.  All  data  have  been  smoothed 
using  a  Gaussian  filter.  The  Gaussian  factor  was  selectable 
and  was  chosen  as  2.5.  The  4096  data  points  taken  at  each 
position  represent  instantaneous  velocities  and  are  used  to 
form  a  probability  density  function  (PDF) .  This  PDF  is 
assumed  to  be  a  Gaussian  distribution.  All  data  outside  2.5 
standard  deviations  of  the  mean  are  eliminated,  and  a  new 
distribution  is  generated.  All  data  lying  outside  2.5 
standard  deviations  for  the  new  PDF  are  eliminated.  A  new 
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figure  6.6  Axial  Velocity  Profiles 
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Figure  6.7  Transverse  Velocity  Profiles 
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Figure  6.10  Covariance  Profiles 


PDF  is  created,  and  the  process  is  complete.  Figures  6.6 
through  6.10  show  how  each  of  the  flow  variables  change  with 
axial  position.  Test  conditions  are  given  in  the  header  of 
each  plot  (refer  to  Table  I) .  The  data  are  presented  here 
to  provide  a  summary  of  results.  Discussion  of  the  results 
will  be  given  as  appropriate  later  in  this  section. 
Attempts  to  perform  spectral  analyses  of  the  data  were 
unsuccessful  because  seeding  rates  were  too  low  at  nearly 
all  conditions.  The  Mie  Scattering  Method  would  have  to  be 
used  to  obtain  the  information  about  vortex  shedding 
frequencies. 

Mie  Scattering 

The  Mie  Scattering  Method  (MSM)  was  used  to  determine 
the  shedding  frequencies  of  the  large-scale  structures.  The 
Mie  scattering  results  are  summarized  in  Figure  6.11. 
Shedding  frequencies  are  shown  for  z  ~  60  mm.  Shedding 

rates  were  the  same  at  all  locations  (for  a  given  U ^ ) , 
except  where  z  =  120  mm  and  uj  =  52.8  cm/s,  where  no 
distinct  frequency  is  observed.  Discussion  of  these  results 
will  follow. 

Discussion 

To  assess  the  character  of  the  bluff  body/ jet 
interaction,  we  first  look  at  the  results  from  the  LDA 
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Figure  6.11  Vortex  Shedding  Frequencies 


measurements.  The  statistical  quantities  obtained  from  the 
LDA  data  may  provide  important  information  on  the 
interaction.  The  mean  velocity  profiles  are  important  since 
they  will  identify,  on  the  average,  the  flow's  behavior. 
The  ability  to  associate  these  mean  flows  with  mean  flows 
from  different,  simpler  geometries  could  provide  useful 
insight  into  the  important  factors  governing  the  flow.  The 
mean  profiles  for  axial  and  transverse  velocity  are  given  in 
Figures  6.6  and  6.7,  respectively.  For  the  axial  velocity 
cases,.  the  mean  flows  above  the  centerbody  at  z  =  6  mm 
resemble  jet  flows.  By  z  *  ^0  mm,  nowever,  the  flows 
(except  for  Uj  =  51.8  cm/s)  resemble  wake  flows.  Sato 
(1960)  suudied  various  slot  jets  (without  anm  lus  flows)  and 
got  mean  profiles  like  those  shown  in  Figure  6.12.  He  also 
observed  that  the  mear  velocity  distrifc  it  ion  changed 
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Figure  6.12  Mean  Velocity  Distribution  (Sato, 1960) 


drastically  around  z/D  of  6  or  7.  Figures  6.6  and  6.7  show 
behavior  very  similar  to  this,  except  these  profiles  are 
affscted  by  the  annulus  flow.  Transverse  scans  at  z  =  10  mm 
were  not  added  on  Figures  6.6  through  6.10  because  the  3-D 
plot  lost  some  of  its  clarity.  The  case  where  z  -  10  mm 
does  very  closely  resemble  the  case  for  z  *  6  mm  for  all 
conditions.  With  the  similarity  seen  in  this  comparison, 
two  definitions  of  flow  "regimes"  will  be  made.  The  "jet 
field"  is  considered  that  region  above  (axial  direction)  the 
bluff  body  up  to  z/D  of  6  or  7.  The  "annulus  field"  is  that 
region  above  the  "jet  field."  It  is  important  to  note  that 
the  "annulus  field"  is  not  just  the  region  above  the 
annulus.  The  "jet  field"  is  the  region  where  the  jet's 
influence  is  strong,  and  the  mean  behavior  of  the  flow  is 
similar  to  that  of  a  jet.  The  "annulus  field"  is  the  region 
where  the  influence  of  the  annulus  is  strong,  and  the  flow 
behaves  like  a  wake.  These  definitions  will  be  useful  in 
later  discussions. 

It  was  hoped  that  the  root  mean  square  of  the 
fluctuating  velocities  and  the  covariance  could  be  used  to 
help  better  define  flow  characteristics.  As  discussed 
earlier  in  the  chapter,  equation  (6.2)  could  be  i'sed  to 
estimate  angular  rotation  rate  for  vortices  assumed  to  be 
solid  body  rotators.  Estimates  of  rotation  rate  for  the 
cases  where  Uj  =  27.8  and  33.3  cm/s  are  given  in  Figure 
6.13.  The  vortex  radii,  r,  were  determined  using 
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Figure  6.13  Estimation  of  Vortex  Angular 
Rotation  Rate 

measurements  based  on  the  photographs,  was  taken  from 

the  data.  The  curves  in  Figure  6.13  provide  a  qualitative 
trend  about  the  flowfiald  which  may  provide  some  insight 
into  the  bluff  body/slot  jet  interaction.  The  covariance, 
shown  in  Figure  6.10,  was  typically  very  small  (absolute 
values  less  than  0.05),  so  this  parameter  was  not  used  to 
help  better  define  the  f.lowfields.  No  trends  or  distinct 
shifts  in  the  parameter  were  seen,  so  little  time  was  spent 
evaluating  the  impact  of  the  covariance. 

As  just  stated,  however,  the  assumptions  used  to  arrive 
at  Eq.  (6.2)  may  be  able  to  help  understand  the 
jet/annulus  interaction.  Their  significance  may  be  brought 
out  in  the  following  way.  The  following  development  has 
been  proposed  by  Hitchcock  (1987)  .  Consider  the  fully 
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developed  duct  flow  shown  in  Figure  6.14.  The  angular 
rotation  rate  of  the  fluid,  given  by  the  curl  of  the 


velocity,  Is 

Wz  " 


1 

2 


Knowing  that 


u 

u 


one  finds  that 


3  y  u 
2  b5 


(6.3) 


(6.4) 


(6.5) 


Since  the  maximum  rotation  speed  will  occur  at  y  ■  b,  Eq. 
(6.5)  can  be  rewritten 


wz  max 

If  one  then  assumes  u 


3  u 
2  b 

U j ,  Eq.  (6.6)  becomes 


(6.6) 


wz  max 


D 


(6.7) 
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Defining 


wz  max 
r  "  ““ 

one  can  see  that  an  estimate  of  angular  rotation  rate  can  be 
made  based  upon  the  "history H  of  the  flow  (i.e.,  before  the 
flow  exhausted  into  the  ambient  stream) .  This  result  is 
interesting  for  several  reasons.  The  values  of  uz  from  Eg. 
(6.7)  are  generally  within  a  factor  of  two  of  those  shown  in 
Figure  6.13.  This  result  is  encouraging  considering  the 
assumptions  made  to  arrive  at  both  results.  (  Whether  the 
flow  is  fully  developed  in  this  investigation  is  a  question 
which  may  be  raised,  but  the  exit  profiles  shown  in  Figure 
6.6  show  mean  velocities  that  tend  toward  a  fully  developed 
duct  flow. )  One  can  also  see  that  there  is  no  dependence  on 
viscosity  in  Eg.  (6.8).  For  a  developing  flow,  viscosity 
would  be  important,  so  parameters  such  as  Reynolds  number 
become  important.  This  is  an  important  point  when 
considering  Strouhal  number,  which  will  be  discussed  next. 

It  was  mentioned  in  Chapter  I  that  use  of  a 
parameter  such  as  St  or  Re  may  be  difficult  in  the 
present  investigation.  To  consider  strouhal  number, 
however,  several  factors  point  to  the  fact  that  Uj  and  slot 
width  D  would  be  the  appropriate  parameters  for 
non-dimensionalizing  the  vortex  shedding  frequencies. 
Even  thoi'gh  the  mean  velocity  profiles  have  shown  a 
combination  of  jet  and  wake-like  flows  (with  the  exception 
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(6.8) 
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of  Uj  -  51.8  cm/s),  the  vortex  shedding  frequency  was  the 
same  in  the  jet  *”d  annulus  "fields."  These  facts 
indicate  that  one  is  justified  in  assuming  that  the  slot 
width  and  jet  exit  velocity  would  be  appropriate  parameters 
with  which  to  non-dimensional ize  the  shedding  frequencies 
(given  in  Figure  6.11).  The  Strouhal  numbers  for  this 
configuration  are  plotted  against  jet  Reynolds  numbers  (Re^) 
in  Figure  6.15.  The  value  of  St  is  nearly  constant  and  has 
the  value  0.098.  The  data  Sato  (1960)  obtained  for  much 
higher  Re  (1500-8000)  produced  a  Strouhal  number  of  0.14. 
This  result  shows,  as  Sato's  did,  that  Strouhal  number  is 
independent  of  jet  Reynolds  number  (for  the  range  of  Re^ 
investigated,  with  the  exception  of  Uj  «  51.8  cm/s).  If  one 
uses  the  definition  of  St,  given  by  Eq.  (1.1),  and  the 
approximation  of  given  by  Eq.  (6.8),  it  is  easy  to  show 
that 


St 


f 


2  f 


(6.9) 


which  says  that  vortex  shedding  frequency  f  is  proportional 
to  the  viscous  rotation  rate  in  fully  developed  duct  flow. 
The  absence  of  Reynolds  number  dependence  again  leads  to  the 
fact  that  the  flow  in  the  duct  is  fully  developed.  It  is 
quite  interesting  that  the  flow  history  is  so  important. 
Recall  that  in  Chapter  I  most  of  the  articles  that  dealt 
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Figure  6.15  Variation  of  Strouhal  Number 
with  Reynolds  Number 

with  the  jet's  near  field  focused  on  the  jet  boundary.  This 
is  certainly  important,  but  the  relation  between  the  flow  in 
the  duct  and  the  flow's  evolution  is  certainly  an 
interesting  one.  Having  analyzed  the  data  obtained  from  the 
LDA  and  the  Mie  scattering  and  having  closely  examined  the 
flow  visualizations,  a  description  of  the  flowfield  dynamics 
can  now  be  presented. 

The  following  description  of  the  flowfield  is 
presented.  For  velocity  ratios  Uj/Uft  up  to  approximately 
1.45,  there  exist  a  "jet  field"  and  an  "annulus  field." 
Beyond  this  velocity  ratio,  the  "jet  field"  extends  into  and 
eliminates  the  "annulus  field."  The  source  of  the 
large-scale  structures  is  the  slot  jet  and  not  the  annulus 
flow.  Evidence  to  support  this  description  will  now  be 
given. 
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Strouhal  number  based  on  the  jet  velocity  and  slot 
width  was  constant  (for  all  cases  except  where  -  51.8 
cm/s) ,  and  this  supports  the  fact  that  the  jet  is 
responsible  for  producing  the  instability  that  results  in 
vortex  formation.  If  the  annulus  flow  was  responsible  for 
vortex  production,  one  would  expect  to  see  a  constant 
Strouhal  number  based  on  annulus  velocity  and  an  appropriate 
length  scale  (such  as  the  width  of  the  bluff  body)  because, 
for  low  Re,  this  flow  should  be  similar  to  shedding  in  the 
wake  of  a  circular  cylinder  (a  flow  known  to  produce 
constant  St  for  low  Re) .  Annulus  Strouhal  numbers  for  this 
experiment,  however,  are  not  constant  since  U  and  D  remain 

a 

constant  while  f  varies.  Also,  the  "stable"  configuration 
for  a  von  Harman  vortex  street  requires  that  b/a  «  0.281, 
where  b  is  the  transverse  separation  between  vortex  centers, 
and  a  is  the  axial  distance  separating  the  centers. 
Evaluation  of  the  quantity  b/a,  based  on  measurements  of  the 
structures  in  the  photographs,  did  not  yield  values  of 
0.281.  Deviations  from  this  were  sometimes  as  high  as  100 
percent,  so  the  inaccuracy  involved  with  measuring  could  be 
ignored . 

The  flow  patterns  seen  in  the  photographs  of  Figure  6.4 
do  show  stable  configuration  of  the  vortices.  The  vertical 
sheet  light  extends  210  mm  in  the  axial  direction,  so  one 
can  view  up  to  70  slot  widths.  This  stability  in  the 
"annulus  field"  is  a  result  of  the  co-flowing  stream's 
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influence  on  the  jet.  Looking  at  Figure  6.4,  the  "trunks" 
of  the  flows  are  like  the  one  depicted  in  Figure  6.16(a). 
Figure  6.16(b)  shows  the  stable  configuration  for  a  plane, 
free  shear  layer  that  was  studied  by  Lasheras  et  al  (1986) . 
It  is  easy  to  imagine  that  a  slight  distortion  (offset  of 
alternating  vortices)  of  the  pattern  in  Figure  6.16(a)  would 
yield  the  configuration  in  (b) ,  which  is  the  stable 
configuration  for  the  2-D  centerbody.  The  geometry  and  flow 
conditions  will  dictate  the  precise  shape  of  the 
configuration.  The  general  shape,  however,  will  remain  the 
same,  and  is  it  the  annulus  flow  that  is  responsible  for 
this  stabilization. 

The  "jet  field"  extends  to  z/D  of  approximately  6  or  7. 
It  is  here  a  significant  redistribution  of  the  mean  axial 
velocity  occurs.  Structures  formed  in  the  "jet  field" 
propagate  through  the  annulus  field  without  much  change. 
However,  as  jet  velocities  increase,  there  is  a  point  where 
the  jet  flow  begins  to  dominate  the  dynamics  of  the  entire 
flowfield.  Looking  at  Figure  6.4  for  =  51.8  cm/s,  one 
may  be  able  to  see  the  interaction  of  some  of  the  vortices 
in  what  should  be  the  "annulus  field."  Near  the  top  of  the 
sheet  light  there  appears  to  be  a  break  in  the  asymmetric 
pattern.  This  is  the  result  of  vortex  interaction  (this  was 
very  clear  when  one  viewed  the  high  speed  movies  of  the 
flow) .  Vortex  interaction  is  not  expected  when  the  annulus 
flow  provides  a  stable  configuration  for  the  vortices,  but 
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the  jet  now  has  sufficiant  energy  to  overcome  a  large 
portion  of  the  annulus'  influence.  The  jet's  strength  has 
now  been  able  to  eliminate,  or  at  least  subdue,  the  "annulus 
field." 

Determination  of  the  approximate  point  where  the 
"annulus  field"  disappears  was  made  possible  by  reanalyzing 
Mie  scattering  data  of  Tankin  (1987)  for  40.4  cm/s  <  < 
71.8  cm/s.  The  slight  discontinuity  in  St  for  -  51.8 
cm/s  given  in  Figure  o.l4  indicated  the  change  probably 
occurred  between  37.0  and  51.8  cm/s,  and  Tankin 's  data 
covered  this  range.  Figures  6.17(a)  through  6.17(e)  show 
the  dominant  frequencies  (at  z  -  80  mm)  in  the  flow.  The 
peak  frequency  at  each  condition  is  given  in  Figure  6.18. 
Vortex  shedding  frequencies  from  this  investigation  are 
repeated  to  provide  a  comparison.  Stepwise  transitions  in 
frequency  indicate  a  jet  flow  is  dominating  the  region. 
This  stepwise  behavior  has  been  reported  for  various  jet 
configurations  by  several  authors,  most  recently  Chen  and 
Schmoll  (1987).  The  frequency  spectra  in  Figure  6.17 
provide  a  beautiful  example  of  the  competing  frequencies  in 
the  jet  over  the  range  of  investigation.  The  dominant 
frequency  remains  dominant  for  relatively  large  jet  velocity 
ranges  before  it  eventually  is  overcome  by  another 
frequency.  It  is  also  interesting  to  note  that  the 
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Figure  6.17  Frequency  Spectra  from  Tankin  (1987) 
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Figure  6.17  Frequency  Spectra  from  Tankin  (1987) 
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Figure  6.17  Frequency  Spectra  from  Tankln  (1987) 


■  Praoant  lnv«»Hg«Mon 
•  From  Tonkin  (18*7) 


#  • 


•  B 


Sam* 

In 

Both 


40  *0  80  100  120  140  180 

Jat  Raynalrfo  Numfcar 


Figure  6.18  Vortex  Shedding  Rate  for  the 
Two-Dimensional  Centerbody 


frequency  which  will  eventually  dominate  the  flow  appears 
well  before  it  becoees  dominant ,  thus  providing  a  hint  of 
what  will  soon  follow  if  jet  velocity  is  increased. 

At  Uj  >37.0  cm/ a  in  the  present  experiment,  the 
shedding  frequency  was  12.125  Hz.  This  condition  was  also 
tested  by  TanXin,  and  tha  same  result  was  obtained.  At  40.4 
cm/m,  a  significant  change  in  frequency  has  occurred. 
Shedding  frequency  is  now  15.125  Hz,  a  value  which  remains 
the  peak  frequency  until  59.2  cm/s.  Because  of  the 
significant  change  at  40.4  cm/m,  it  was  assumed  that  this 
signified  the  initiation  of  the  jet's  dominance  over  the 
entire  flovfield.  The  approximate  velocity  ratio  Uj/Ua  at 
which  this  dominance  is  seen  is  1.45.  Velocity  ratios  less 
than  this  should  produce  jet  and  annulus  "fields,"  while 
ratios  greater  than  1.45  should  result  in  jet  domination. 
This  description  does  not  ignore  the  influence  of  the 
annulus  flow;  its  main  point  is  simply  that  the  "annulus 
field"  has  disappeared  and  lost  some  measure  of  its 
stabilizing  strength. 

Having  completed  the  description  of  the  flowfield 
dynamics,  how  ccn  it  be  used  to  explain  the  observed 
bifurcations?  Viewing  Figure  6.4,  one  can  see  some  patterns 
which  are  rather  easily  explained.  At  Uj  *  18.5  cm/s,  we 
see  the  strong  influence  of  the  annulus  flow  as  it  tries  to 
wrap  jet  fluid  around  itself,  thus  forming  the  large  vortex 
structures  which  are  oriented  "up."  At  -  51.8  cm/s,  the 
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j«t  fluid  wraps  ambient  fluid  around  itsalf  and  vorticaa 
with  a  "down"  orisntation  art  foraed.  Tha  flow  pattam 
change  at  -  22.2  ca/s  rssults  bscauss  of  tha  intaraction 
of  tha  annulus  and  jat  flows  whan  thair  mean  axit  valocitias 
ara  aqual.  Pigura  €.19  shows  that  thara  is  a  diffaranca 
batwaan  thaoracical  and  actual  aaan  valocitias.  At  thasa 
low  flow  rates ,  it  is  axpactad  that  tha  flow  profile  would 
daviata  fro*  tha  idealized  (1-D)  casa.  At  -  22.2  cm/ a, 
aaan  axit  valocitias  for  tha  annulus  and  jat  ara  aqual,  and 
this  results  in  tha  strange  flow  pattern.  Tha  transverse 
spread  of  tha  jat  is  vary  saall.  It  appears  like  tha  jat  is 


Figure  6.19  Axial  Velocity  Profile,  Transverse  Scan 
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confined  to  that  apace  directly  above  the  bluff  body.  The 
TiOj  particles  sees  to  act  almost  like  a  solid  wall  meant 
to  separate  the  two  annulus  flows.  With  equal  velocities, 
little  entrainment  is  expected,  so  the  transverse  spread  of 
the  jet  is  minimised. 

The  bifurcation  at  -  37.0  cm/s  is  quite  interesting. 
This  condition  is  very  close  to  the  point  where  the  "annulus 
field"  disappears.  The  nearly  straight  line  pattern  from 
Tankin  (1987),  shown  in  Figure  4.1(h),  la  not  Batched 
exactly  in  the  present  experiment  (see  Figure  6.4(e)),  but 
the  bifurcation  can  still  be  described.  At  this  jet 
velocity,  scan  uaall  disturbances  are  seen  in  the  form  of 
small  white  "dots"  which  should  be  the  beginnings  of 
largo-scale  structures.  These  disturbences,  however,  are 
very  weakly  amplified.  Inspection  of  the  mean  axial 
velocity  profile  at  s  -  60  mm.  Figure  6.20,  will  help 
explain  why  this  occurs.  The  flat  velocity  profile  says 
there  is  no  shear  in  the  "annulus  field."  Therefore,  the 
disturbance  (which  results  from  the  jet  flow  but  whose 
evolution  is  still  governed  by  the  annulus  flow)  is  not 
aaplified.  If  jet  velocity  is  increased  slightly  (so  U^/Ua 
>  1.45),  amplification  of  the  disturbance  will  occur,  and 
the  jet  flow  will  dominete  the  entire  flowfield. 

The  flow  pattern  change  at  Uj  -  22.2  cm/s  is  important, 
but  the  bifurcation  at  ■  37.0  cm/s  is  critical.  While 
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Figure  6.20  Axial  Vatoctty  Profile,  Tranavarea  Scan 
z  -  60  mm 

the  former  bifurcation  point  occur*  bacausa  of  velocity 
counterposition,  tha  lattar  represents  a  point  around  which 
tha  governing  factors  of  tha  flovf ield  will  change.  Up  to 
Uj  ■  37.0  cm/s,  one  expects  tha  structures  formed  in  tha 
Njet  field"  to  be  carried  by  the  stabilising  annulus  flow 
without  much  change.  Beyond  this  velocity,  however,  jet 
strength  is  sufficient  to  subdue  the  influence  of  the 
annulus  flow,  and  a  jet-dominated  flow  results.  Proa  the 
data  taken,  the  flow  pattern  change  at  22.2  cm/s  did  not 
reveal  a  change  in  the  character  of  the  flow  similar  to  the 
one  seen  at  37.0  cm/s. 
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There  is  on*  more  point  which  can  be  »ddr*eeed.  Flow 
patterns  are  always  asymmetric  (i.e.,  position  of  ths 
struct urss  sltsmatss  fro*  sids  to  sids) .  Two  explanations 
ars  offsrsd.  First,  ths  analogy  of  an  axisymmetric  jat 
which  is  being  squeexod  nay  ba  used.  Ths  flow  patterns 
during  a  simulated  contraction  process  are  shown  in  Figure 
6.21.  The  vortices  in  (a)  start  to  get  closer  and  closer  as 
the  jet  width  narrows,  At  a  certain  point,  (c) ,  the 
vortices  can  no  longer  be  brought  closer,  since  the 
vortices  do  not  interact  (coalesce) ,  they  oust  align 
theaselves  as  in  (d)  .  The  geometry,  then,  would  specify 
whether  or  not  symmetric  patterns  would  be  observed. 
Changing  geometries  (a.g. ,  slot  or  bluff  body  widths)  could 
change  the  patterns.  This  simplified  analysis  does  not 
account  for  end  effects  which  will  play  an  important  role  in 
determining  if  the  vortices  can  be  shed  symmetrical ly . 

The  second  explanation  involves  a  feedback  mechanism. 
According  to  the  description  of  the  flovfield  dynamics,  the 
configuration  shown  in  Figure  6.16(a)  is  stable.  This 
asymmetric  pattern  will  set  up  an  induced  pressure  field 
that,  even  though  it  is  very  weak,  can  affect  the  flow  at 
the  jet  exit.  Ho  and  Nosseir  (1962)  used  a  jet  impinging  on 
a  wail  to  show  that  resonance  was  established  based  on  a 
feedback  mechanism.  A  symmetric  vortex  shedding  is  unstable 
in  the  "annulus  field,"  so  asymmetric  patterns  are  formed 
which  create  the  induced  pressure  field  that  results  in 
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Figure  6.21  Conception  of  the  Reason  for  Asymmetric 
Vortex  Shedding  In  a  Thin  Slot  Jet 


ityiMtric  shadding  at  tha  jat  axit.  strong  avidnnca  for 
tha  axistanca  of  this  aachanisa  in  this  flow  was  aaon  in  tha 
Mia  aoattaring  axpariaanta.  For  ■  33,3  cm/%,  tha 
aaasuraaant  apot  was  looatad  in  a  rag ion  (>  ■  20  wm)  whara 
no  structuras  ax is tad,  although  ona  could  aaa  that  tha 
" intar faca  linaa"  vara  vavaring.  Tha  fraquancy  of  this 
vavaring,  aa  datarainad  by  Mia  acattaring  and  shown  in 
Figura  6.22,  was  11.375  Hs.  This  is  tha  saaa  fraquancy 
datarainad  whan  tha  data  vara  takan  in  a  ragion  vhara 
structuras  did  axiat  (a  ■  40  aa  and  s  -  60  aa) .  This 
providas  strong  avidanca  for  tha  axistanca  of  a  faadback 
aachanisa  in  tha  flows  undar  invastigation. 
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VII*  Concilia  Iona  and  BlfiflUDflatlflM 


Tha  character  of  tha  interact ion  batwaan  bluff  body  and 
jat  f Iowa  haa  baan  datarainad.  Tha  concapt  of  jat  and 
annul ua  "fialda"  can  ba  uaad  to  axplain  tha  intaractiona  up 
to  valocity  ratioa  U^/U^  of  1.45.  Bayond  this  valua,  tha 
flowfiald  dynaaica  ara  govarnad  by  tha  jat,  but  tha  annulus 
flow  cannot  ba  ignorad.  Tha  jat  datarainoa  tha  shadding 
fraquancias  of  tha  largo-ccala  atructuras  and  govarno  tha 
flow  in  tha  ragion  abova  tha  bluff  body  to  s/D  of 
approxiaataly  6  or  7  (i.a.,  tha  "jat  fiald").  Tha  "annulus 
fiald"  is  that  ragion  abova  tha  "jat  fiald"  wharain  the 
annulus  flow  craatas  a  stabla  vortax  configuration  and 
govams  tha  flow's  dynaaica. 

Tha  two  obaervad  flow  pattarn  changas,  bifurcations, 
vara  partially  axplainad.  Tha  flow  pattarn  at  IK  -22.2 
ca/s  rasultad  froa  flow  conditions  whara  jat  and  annulus 
valocitias  wars  naarly  idantical.  Tha  bifurcation  at  - 
37.0  ca/a  was  a  critical  one  sinca  it  idantifiad  tha  point 
abova  which  tha  "annulus  fiald"  haa  baan  aliainatad.  Thara 
ara  aany  iaportant  factors  in  tha  intaraction  of  tha  jat  and 
annulus  flows,  and  flow  pattarn  changas  hava  baan  partially 
axplainad  by  using  various  labals.  Thara  ara  still  aany 
datails  of  tha  flowfiald  which  could  ba  uncovarad. 

Racoaaandations  for  furthar  rasaarch  would  ba  to 
invaatigata  tha  intaraction  at  highar  valocitias  (ona  ordar 
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of  aagnitude  at  least)  to  determine  if  the  same  daacription 
ia  valid.  Redesign  of  tha  cantarbody  ahould  firat  ba 
accoaplished,  however.  Also,  testing  ahould  ba  dona  in  a 
low  turbulanca  facility  to  avoid  tha  problems  with  the 
annulua  turbulanca  lavala. 

No  point  vhara  tha  annulua  flow  dominates  tha  antira 
flowfiald  waa  obaerved  in  thia  invastigation.  The 
visualizations  froa  Tankin  (1987),  hovavar,  aay  auggaat  that 
thara  ia  a  point  vhara  thia  doaa  occur.  Figura  4.1(a)  shows 
tha  condition  vhara  -  4  cm/m.  No  structuras  ara  aaan 
aaanating  froa  tha  slot.  As  jat  valocity  is  incraasad, 
structuras  foraad  at  tha  slot  exit  bacoaa  evident.  The  flow 
at  Uj  ■  4  ca/s  is  auch  like  tha  flow  behind  a  bluff  body  (so 
the  slot  jat  only  acts  as  an  orifice  to  provide  tha  seed 
needed  to  visualize  tha  flow) .  Expariaants  at  vary  low 
Uj/U^  night  provide  information  concerning  the  interactions 
of  the  annulus  and  jat  flows  where  tha  annulus  doainatas  tha 
antira  flowfiald. 

Thara  ara  aany  othar  araas  vhara  thesa  rasults  could  ba 
applied.  In  Computat.u  nal  Fluid  Dynaaics,  these 
incoapressible,  2-D  flows  could  ba  aodelled,  and 
those  rasults  could  ba  coaparad  with  those  froa  this 
axpariaant.  For  this  reason  all  data  collected  during  tha 
LDA  aaasuraaants  ara  presented  in  Appendix  C.  Availability 
of  data  for  a  unique  configuration  is  iaportant,  so  these 
data  ara  included  in  this  thesis.  Critical  point  theory 
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could  be  applied  to  further  analyze  these  flovfieldo. 
linear  and  non-linear  instability  analyses  night 
provide  insight  in  flow  regines  where  tractable. 
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APPENDIX  A:  DETERMINATION  OF 
tici4  CONCENTRATION 

Flow  visualization  used  in  this  sxpsrissnt  has  been 
described  in  detail  (Chan  and  Roquemore ,  1986) .  Titanium 
tatrachlorida  (Ticl^)  is  introducad  to  a  dry  fluid  and 
axhausts  as  a  sixtura  through  tha  axparisantal  apparatus 
( jet,  ate) .  Whan  tha  TiCl4  contacts  vatar  vapor  (H20) , 
titanium  dioxida  (Ti02)  is  forsad.  This  raaction  occurs 
undar  naarly  isotharaal  conditions,  so  ralativaly  little 
energy  is  lost  or  gained  by  tha  flow.  Tha  subnicron  TiO^ 
particles  respond  to  flow  variations  naarly  instantaneously. 
Reactive  flow  visualization  techniques  like  tha  one  described 
above  provide  important  inuight  to  tha  processes  occurring  in 
various  flows. 

An  important  question  concerning  this  flow  visualization 
technique  had  not  been  answered:  NNhat  is  the  seed 
concentration?"  Density,  and  acre  importantly  density 
gradients,  are  critical  factors  in  jet  and  wake  flows; 
therefore,  the  seed  concentration  must  be  small  so  that 
density  gradients  are  not  introduced.  A  model  to  determine 
seed  concentration  for  flow  visualization  has  been  developed. 

A  typical  experimental  set-up  is  shown  in  Figure  A-l. 
Dry  air  passes  through  a  flow  meter  and  then  to  a  seeding 
vessel.  The  seeding  vessel  is  shown  in  greater  detail  in 
Figure  A-2.  Flow  enters  and  can  be  bypassed  to  any  degree 
desired  (0-100  percent) .  Flow  enters  the  container,  wherein 
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liquid  TiCl4  is  stored,  through  three  1/4  inch  lines.  Each 
of  these  lines  has  a  ball  valve  so  that  flow  say  go  through 
1,  2,  or  3  of  the  lines.  Many  different  seeders  have  been 
built;  most  give  adequate  seeding  rate  with  only  one  input 
line.  This  is  what  was  used  in  this  analysis.  The  lines 
inside  the  container  are  angled  at  -45  degrees  fros  the 
horizontal  and  are  canted  so  the  flow  leaves  nearly  parallel 
to  the  wall.  The  design  enhances  mixing  and  ensures  adequate 
amounts  of  TiCl4  get  into  the  flow. 

Mass  transfer  theory  can  be  used  to  determine  the  seed 
concentration.  The  seeder  is  designed  to  get  the  liquid 
TiCl4  rotating.  This  has  been  observed  in  a  clear  plexiglass 
seeder  designed  for  water  vapor  uptake.  Mass  transfer  for  a 
laminar  rotating  disk  has  been  analyzed  (Levich,  1962)  and 
was  used  to  model  the  flow  in  the  seeder.  The  mass  flux  from 
the  surface  to  the  fluid  is  given  by 

j  -  k  (cs  -  cQ)  -  .62  Dftb2/3  j/-V6wl/2  ^  (A.  1) 

where  c  is  the  surface  concentration  of  TiCl.,  and  c  is  the 
overall  concentration  of  TiCl4.  U)  is  the  angular  rate  at 
which  the  disk  rotates  (or  the  angular  rate  at  which  the  flow 
passes  the  liquid  surface) .  D&b  is  the  binary  diffusion 
coefficient.  The  value  of  V ,  the  kinematic  viscosity,  is 
taken  as  that  for  air.  If  we  assume  cQ/cm  «  1,  we  get 

k  -  .62  Dafa2/3  y“1/6GJ1/2  (A.  2) 
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Proa  o  mm  balance,  wo  also  know  that 


cD  v  -  A  Jt  (es  -  c0) 


(A. 3) 


whoro  V  is  th«  volunotric  flow  rato  of  air,  and  A  Is  the 
surface  area  of  the  liquid  TiCl4.  The  concentration  of 
TiCl4,  in  parts  per  thousand  by  voIum,  is  thus 


c  A  k 

-2  .  ( - )  1000  (A. 4) 

C  V  +  Ak 


In  order  to  determine  surface  concentration  we  use 


(A.5) 


where  is  the  vapor  pressure  of  TiCl4,  and  R  is  the 
universal  gas  constant.  Several  methods  for  determining 
vapor  pressure  are  available  and  are  summarised  by  Reid  et  al 
(1977).  One  method  was  chosen  that  was  simplest  for  computer 
application,  although  it  was  not  designed  as  such  (Reidel, 


1954) . 


ln  pvp  r  "  A+  "  <B+/Tr)  +  C+  ln  Tr  +  D+  Tr6  (A’6) 

where  the  subscript  r  denotes  a  reduced  property  ( i . e . ,  one 
normalized  by  its  appropriate  critical  property) .  We  have 
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-  -35  Q 

(A. 7a) 

-  -36  Q 

(A.  7b) 

-  42  q  +  a. 

(A. 7c) 

-  -o  c 

(A. 7d) 

-  . 0838  (3.758  -  Ot  „) 

c 

(A.8) 

0.315  \ph  +  In  Pc 

i 

(A.9) 

0.0*3*  1 pb  -  In  Tbr 

-3S  +  (3«/Tbr)  +  43  in  Tbr  -  Tbr 

(A.  10) 

In  order  to  determine  the  vapor  pressure  using  equation 
(A. 6) ,  critical  properties  for  TiCl4  must  be  known.  Boiling 
temperature  was  also  required  and  was  obtained  from  the 

Ohficfli  -Bnalniiri* _ HanflbggX  (1973).  Critical  temperature 

and  pressure  of  TiCl4  were  determined  using  the  method  of 
Lydersen  (1955).  Critical  properties  are  determined  using 
structural  contributions.  We  have 


Tc  "  Tb  <0*567  +  f<At>i*  (A*11) 
Pc  -  M  [0.34  +  <£(Ap)i“2]  (A.  12) 

where  Tc  is  in  degrees  Kelvin,  and  Pc  is  in  atmospheres.  M 
is  the  molecular  weight,  and  At  and  Ap  are  the  structural 
contributions.  These  values  are  given  in  Table  A-I  for  Si 
and  Cl.  No  data  could  be  found  for  TiCl4  or  Ti,  so 
approximations  were  made.  Due  to  bonding  similarity,  use  of 
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Si  viluti  should  yield  acceptable  results.  The  desired 
answer  need  not  be  perfect,  so  little  tine  was  spent 
searching  for  a  better  aodel. 


TABLE  A- I  Lydersen's  Critical  Contributions 
(Lydersen,  1955) 


Increments 

*t 

AP 

—  Cl 

0.017 

0.320 

1 

—  Si  — 

0.03 

0.54 

1 

Using  equations  (A. 11)  and  (A. 12),  one  can  find  the  critical 
properties  for  TiCl4.  The  important  properties  of  Tici4  are 
listed  in  Table  A-II.  Equations  (A. 6)  through  (A. 10)  are  now 
used  to  determine  vapor  pressure  after  a  temperature  has  been 
specified. 


TABLE  A-II  TiCl4  PROPERTIES 


N  (91/sole) 
Tw  <fK) 


3 

P=  <«tB> 


<°1C) 


189.71 

409 

624 

40.5 


Data  for  vapor  pressure  of  TiCl4  existed  in  the  Chemical 
Engineers 1  Handbook  (1973).  These  data  and  the  model  using 
(A. 6)  through  (A. 10)  are  shown  in  Figure  A-3.  A  correction 
to  the  model  was  made  to  better  natch  the  data  from  the 
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Chemical _ Enalnaara 1  Handbook .  Sine*  vapor  praaaura  nodal a 

•m  conplax,  it  was  dacidadly  aaaiar  to  usa  tha  Raidal  nathod 
and  corract  it  Instaad  of  trying  to  curva  fit  tha  axisting 
data.  Tha  nodal,  with  ita  corract ion,  ia 


Pyp  *  [l.«5  ♦  0.0046857  (260  -  T)  ]  P^  c  (A.  13) 

which  ia  valid  from  260°K  to  330°K.  „  ia  tha  vapor 
praaaura  calculatad  with  Raidal 'a  nathod.  Tha  graataat  arror 
occurs  at  260°K  and  ia  approxinataly  4.5  parcant.  Tha  arror 
quickly  dacraaaaa  with  taaparatura  incraaaa  and  ranaina 
naarly  constant  at  1  parcant  ovar  tha  raaaindar  of  tha  val.** 
ranga. 
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The  next  step  in  dtttraining  concentration  was  to 
evaluate  tha  diffusion  coefficient,  D^.  Using  the  sethod  of 
Bird  et  al  (1900) ,  va  have 

(f 3 [ ( l/n\ ♦ ( l/Nw) ] > 1/2  cs-cs 

D  -  .0019593 - #— -  ( - }  (A.  14) 

vhere  T  is  t— psrature  in  °K,  P  is  the  pressure  in 
ataoepheres,  and  are  Molecular  weights  of  the  two 

Molecules  (TiCl4  and  air)  in  grans  per  sole,  is  tha 

characteristic  dinension  of  the  Molecules  in  Angst roes,  and 
0  oab  A*  collision  cross  section  for  diffusion 

(diMensionless) .  To  calculate  we  Must  use 


-  o.s  ( am  +  crh) 


(A. 15) 


where  values  of  (J  are  given  for  various  gases  by  Svahla 
(1902).  For  TiCl4,  the  value  of  G  was  aasunad  to  be  that  for 
Sld4  since  no  data  were  available  for  TiCl4.  To  calculate 
0 ,  an  enpirical  relation  for  nonpolar  gases  given  by  Neufeld 
at  al  (1972)  was  used. 


C  E  6 

exp(DT*)  exp (FT*)  exp(HT*) 


vhere 


T* 


k  T 


(A. 16) 


(A. 17) 
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and 

A  -  1.06036 
B  -  0.15610 
C  -  0.10300 
D  -  0.47635 
K  -  1.03567 
F  -  1.52996 
O  -  1.76474 
H  -  3.89411 


k  ir  the  well-known  Doltsmann  constant ,  and  6^  la  an 
intermolecular  force  constant  vhich  is  calculated  as 


««b  -  ( e«  it,)1'*  (a. is) 


Values  for  6  are  given  by  Svehla  (1962) .  Again,  the  value  of 

6  for  TiCl.  was  assumed  the  sane  as  that  for  SiCl. . 

4  4 

In  order  to  calculate  the  mass  transfer  coefficient  in 
equation  (A. 2).  ve  aust  know  the  kinematic  viscosity  of  the 
aixture.  Assuming  the  concentration  of  TiCl4  is  low 
everywhere  except  very  near  the  liquid,  ve  need  only 
calculate  the  kinematic  viscosity  of  air.  We  know 


V  -  ,U/,Q  (A.  19) 


where 


P 

0  - -  (A. 20) 

R  T 

To  calculate  the  viscosity,  we  aust  use  the  equation  from 
Bird  et  al  (1960) 
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(A. 21) 


-  (*T)1/a  qm 

LL  -  2. *€93X10  *  -  < - ) 

*  «-• 

vbtrt  N  is  th«  aoltoultr  weight  of  air  in  grams  per  nolo, 
CT  io  the  characteristic  dimension  of  the  air  aolecule  in 
Angotro— ,  and  Q,#  io  the  collision  cross  section  for 
viscosity.  The  collision  cross  section  is  s  slowly  varying 
function  (near  unity)  which  indicates  how  close  the  soleculea 
under  consideration  behave  like  perfectly  elastic  spheres. 
To  evaluate  ,  we  use  an  enpirical  relation  proposed  by 
Meufeld  at  al  (1972). 


AC  K 

fi/,  “  “S  ♦ - ♦ -  (A.22) 

M  exp (DT*)  exp (FT*) 

where 


k  T 

T*  - -  (A. 23) 

€ 

and 

A  -  1.16145 
B  -  0.14874 
C  -  0.52487 
D  -  0.77320 
E  -  2.16178 
P  -  2.43787 

(T  and  €  can  be  obtained  from  existing  data  for  air  (Svehla, 
1962) . 

We  say  now  calculate  k,  the  mass  transfer  coefficient, 
if  wo  can  evaluate  the  angular  speed  CJL  An  estisate  of  CJ  was 
used. 
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U)  -  0.25  (V/C) 


(*•24) 


vhtn  v  is  tlM  velocity  of  air  loading  tha  1/4  inch  lina 
which  antara  tha  top  of  tha  caniatar  (aaa  Figure  A-2) ,  and  C 
la  tha  oaniatar'a  ins ids  circuafaranea  (TTd). 

Ona  final  adjustment  oust  ba  nada  bacausa  tha  nasa 
transfer  nodal  developed  by  Levich  was  valid  only  for  laminar 
flow.  To  correct  for  turbulence,  a  flat  plate  analogy  was 
used.  Welty  et  al  (1974)  have  given,  for  a  flat  plate. 


klam  “  0.332  (D^x)  Rex1/2  Sc1/3  (A. 25) 

and 

*turb  "  °*02*2  <*>*i/*>  r*x4/5  ScV3  (a-26> 


Forming  a  ratio,  wa  have 


— -£-  -  0.08795  (Raw)0’3  (A. 27) 

v  * 

*lam 

This  factor  varies  from  unity  at  Rex  of  3,300  to 
approximately  5.5  at  Rex  of  1,000,000.  Assuming  tha  particle 
travels  a  distance  d,  tha  diameter  of  tha  canister,  we  find 

V  d 

Re  -  — - -  (A.  28) 

^air 
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wb sra  is  tha  valocity  in  tha  can i star  and  is  assumed  to  bs 
25  psrosnt  of  ths  vslooity  1 saving  ths  1/4  inch  inlst  lins 
(i.s. ,  0.2S  v) . 

Ths  squations  in  this  appsndix  hava  baan  put  in  a 
coaputar  program,  CONCENTRATION ,  in  ordar  to  quickly  analyzs 
ths  TiCl4  concantration  prasant  in  diffarant  flows  undsr 
various  conditions.  Ona  inputs  teaparatura ,  pressure,  and 
jat  valocity  (this  spacifias  voluastric  flow  rata  for  ths 
spacific  tast  arrangaaant) .  Sons  rasults  ara  given  in  Tabls 
A-III.  Tha  concantration  of  TiCl^,  givan  in  parts  par 
thousand  by  volume,  is  saall  in  all  instancas  and  can  bs 
justifiably  ignorad  as  an  influanca  on  tha  flow  undsr  study. 
This  ganaral  stataaant  is  not  applicabla  in  flows  of  vary  low 
valocity  (  <5  cm/m)  and  high  saading.  Changas  in  flow 
pattams  hava  baan  obsarvsd  whan  saading  ratas  incraasad  to  a 
sax i sum  valua  at  thasa  low  vslocitias.  Ths  dsnsity  and 
tasparaturs  gradients  inducad  by  adding  tha  saad  have  a 
longar  tins  ovar  which  to  act,  so  natural  convsction  can  no 
longar  bs  ignorad.  In  othsr  instancas,  howsvsr,  ths  sffset 
of  ths  saad  on  ths  flow  dsvalopssnt  and  transport  of  mass, 
aoaantua,  and  anargy  can  ba  ignorad. 


TABU  A-III  CONCENTRATION  OF  TiCl4 
AXXSYMMETRIC  ual  (D  ■  1  cm) 


TEMP  (°K) 

P  (*tm) 

k  (cm/*) 

C  (ppt) 

260 

1 

7.52 

0.18 

270 

1 

7.52 

0.38 

280 

1 

7.52 

0.76 

290 

1 

7.52 

1.43 

300 

1 

7.52 

2.57 

310 

1 

7.52 

1.25 

4.39 

320 

1 

7.52 

1.27 

7.19 

330 

1 

7.52 

1.28 

11.36 

260 

1 

3.05 

0.55 

0.21 

270 

1 

3.05 

0.56 

0.45 

280 

1 

3.05 

0.57 

0.89 

290 

1 

3.05 

0.58 

1.67 

300 

1 

3.05 

0.59 

3.00 

310 

1 

3.05 

0.60 

5.12 

320 

1 

3.05 

0.61 

8.39 

330 

1 

3.05 

0.62 
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APPENDIX  B:  VALIDATION  OF  THE 
MIE  SCATTERING  METHOD 


A  set  of  experiments  was  designed  to  test  the  validity 
of  the  Mie  scattering  method  (KSM)  for  determining  the 
shedding  frequency  of  vortex  structures  in  various  flows. 
The  peak  frequency  of  velocity  fluctuations  corresponds  to 
the  rate  at  which  vortex  structures  pass  and  can  provide 
important  information  for  understanding  flow  dynamics.  For 
example,  vortex  pairing  has  occurred  when  the  frequency  of 
the  fluctuating  velocity  is  halved  (Chen  and  Schmoll,  1987). 
Frequency  spectra  can  be  obtained  from  LDA,  hot  wire,  or 
hot/ cold  wire  experiments.  When  LDA  is  used,  seeding  rates 
must  be  sufficient  to  ensure  that  data  rates  fulfill  the 
Nyquist  criterion  (i.e.,  the  data  rate  must  be  at  least 
twice  the  highest  frequency  component  sne  is  trying  to 
measure) .  This  must  be  done  to  avoid  losing  any  frequency 
higher  than  the  Nyquist  frequency  (half  the  sampling  rate) 
and  avoid  having  the  undersampled  high  frequencies  show  up 
as  spurious  low  frequencies,  a  phenomenon  known  as  aliasing 
( Malms tad t  et  al,  1981).  As  discussed  in  Section  4,  seeding 
problems  arise  if  flow  rates  are  extremely  low,  such  as 
those  in  the  present  study.  Data  rates  are  often  too  low  to 
provide  spectral  information.  The  MSM  could  provide  a 
viable  way  of  determining  important  shedding  frequencies. 
The  method  was  never  previously  validated,  so  a  set  of 
experiments  was  designed  to  check  its  validity.  Previous 
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us*  of  th*  MSH  for  determining  shedding  fr*qu*nci*s  was 
inconclusive  because  aany  fr*qu*ncl*s  of tan  resulted  from 
the  Fourier  analyses.  Coeplex  signals  made  it  impossible  to 
interpret  the  results  in  a  satisfactory  manner.  If  one 
could  fix  the  vortex  shedding  frequency  at  a  given  value, 
the  MSM  could  be  used,  and  the  results  should  be  much  easier 
to  interpret. 

An  axisymaetric  jet  configuration  was  used  to  test  the 
Mie  scattering  method.  The  test  set-up  is  shown  in  Figures 
B-l  and  B-2 .  A  facility  for  diagnostic  testing  was  being 
established  by  the  project  sponsor,  and  this  is  where  the 
validation  took  place.  The  flow  system  used  dry  air  brought 
through  a  flow  meter  (Fisher  Porter  model  using  Tube  No. 
85-27-10/77) ,  through  a  seeder,  and  then  through  a  1  meter 
tube.  Nozzle  exit  diameter  was  1  centimeter.  Acoustic 
forcing  was  used  in  this  experiment  in  order  to  "drive"  the 
flow.  A  sinusoidal  input,  one  whose  frequency  and  amplitude 
could  be  varied,  is  provided  to  the  speaker.  The  acoustic 
wave  is  transmitted  to  the  flow  through  a  thin  rubber 
membrane  (part  of  a  surgeons  glove)  which  is  placed  over  the 
line  that  fits  into  the  bottom  of  th*  tee  which  it.  located 
at  the  bottom  of  the  1  meter  tube.  Flexible  tubing  was 
inserted  in  the  inner  horn  of  the  speaker  and  led  to  the 
bottom  of  the  tee  fitting.  The  glove  is  used  to  protect  the 
speaker  and  other  equipment  from  potential  contamination 
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Fiber 

Optic 

Cable  Experimental  Table  (Top  View) 
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Photomultiplier 

Tube 


Figure  B-2  Optical  Arrangement  for  Mis  Scattering 
Method  Validation 
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from  the  HC1  which  is  formed  in  ths  rsaction  of  ths  TiCl^. 
Ths  optical  arrangamant  is  tha  sama  as  tha  ona  dascribad  in 
Saction  2. 

A  variabla  spaad  stroba  was  usad  to  "freeze"  tha  motion 
of  tha  jat.  This  visualization  tachniqua  halps  idantify  the 
structures '  shadding  fraquancias.  A  stroba  and  tha  MSM  vara 
usad  to  datarmina  tha  natural  (i.e. ,  unf oread)  shadding 
fraquancy  of  tha  jat  shown  in  Figura  B-l  wi  i  velocity  4 
m/s.  Tha  facility  usad  did  not  allow  for  tha  elimination  of 
outside  disturbances,  so  tha  presence  of  several  different 
frequencies  was  expected.  It  was  difficult  to  exactly 
"freeze"  tha  jat  motion  using  tha  stroba,  and  tha  Fourier 
transform  of  tha  light  scattering  signal,  shown  in  figure 
B-3 ,  did  contain  several  fraquancias.  Tha  flow  appeared 
closest  to  being  frozen  by  tha  stroba  at  approximately  365 
Hz.  Tha  value  of  365  Hz  is  not  seen  exactly  in  Figura  B-3, 
but  tha  two  strongest  peaks  surround  that  value.  Tha  peak 
at  437 , 5  Hz  may  result  from  external  disturbances  or  the  jat 
may  be  nearing  a  flow  condition  where  transition  to  a  higher 
shadding  rata  will  taka  place.  Tha  fraquancy  is  certainly 
not  negligible,  but  tha  MSM  cannot  tall  tha  source  of  tha 
various  instabilities.  This  must  be  left  to  tha 
investigator. 

Acoustic  forcing  was  than  usad  to  raise  tha  level  of 
tha  natural  instability,  assumed  to  be  365  Hz,  wall  above 
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Axisymmetric  Jet 
Velocity  •  4  m/s 
Diameter  -  1  cm 


ill 


A 


950 


Figure  B-3.  Natural  Jet  Vortex  Shedding  Frequency  Spectrum 


Axisymmetric  Jet 
Velocity  -  4  m/s 
Diameter  ■  1  cm 

*  forcing  * 


Figure  B-4.  Forced  Jet  Frequency  Spectrum 


the  background  noise  and  external  disturbances .  This  is 
exactly  what  Crow  and  Champagne  (1971)  did  in  their  classic 
experiment.  If  the  MSM  is  unambiguous ,  only  one  frequency 
should  be  evident.  If  not,  perhaps  the  additional 
frequencies  would  become  clear.  For  example,  one  might 
speculate  that  both  vortex  sise  and  spacing  would  yield 
frequencies  when  spectral  analysis  was  performed.  If  two 
frequencies  result  when  the  jet  is  forced,  one  could  measure 
vortex  sise  and  spacing  and,  with  the  velocity,  see  how 
close  this  is  to  the  observed  result. 

Figure  B-4  shows  the  frequency  spectrum  when  the  jet  is 
forced  at  its  natural  frequency.  Only  one  component  of 
frequency  exists,  and  it  is  the  shedding  frequency  of  the 
large  scale  structures.  This  fact  was  verified  using  the 
strobe  (which  was  synchronised  with  the  oscillator) .  Note 
the  scales  on  the  y-axes  of  Figures  B-3  and  B-4.  Although 
the  vertical  scales  in  the  two  figures  display  arbitrary 
units,  the  scales  may  be  compared  against  one  another.  The 
natural  jet  spectrum  must  be  multiplied  by  a  factor  of  8  to 
achieve  the  levels  shown,  while  the  forced  jet  has  a  factor 
of  only  2.  Thus,  forcing  has  effectively  eliminated  the 
background.  These  experiments  prove  the  MSN  is  a  useful 
tool  in  obtaining  the  shedding  frequencies  for  the  vortex 
structures  in  different  flowfields. 
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APPENDIX  C 


X.  RUN  SUMMARY  OP  THE  6  FLOW  CONDITIONS  IN  THE  PRESENT 
STUDY 

2.  DATA  FOR  ALL  RUNS 

-  SET  UP  AS  2  ROWS  OF  INFORMATION  FOR  EACH  DATA  POINT 

-  SET  UP  IS  AS  FOLLOWS 


(u'v')n 


SU  *U  SV  *v  0  0  0  (T.K.E.  )n 


u  -  Mean  axial  valocity 
v  -  Maan  transvarsa  (radial)  valocity 
u  g  ~  Root  nean  square  of  fluctuating  axial  valocity 
u  s  -  Root  naan  square  of  fluctuating  transverse 
valocity 

(u,v,)n  -  Normalized  Reynold  stresses 
Su  -  Skewness  of  u 

-  Kurt os is  of  u 
Sv  -  Skewness  of  v 

-  Kurt os is  of  v 

(T.K.E. )  -  Normalized  turbulent  kinetic  energy 


Table  C-I  LDA  Run  Sunary 


Ua  (ca/s) 

(cm/ a) 

*  (*■) 

RUM  SEQUENCE 

DATA  QUALITY 

27.5 

18.5 

6 

9029-9054 

GOOD 

10 

9055-9080 

GOOD 

20 

9081-9106 

GOOD 

40 

9107-9132 

GOOD 

60 

9133-9158 

GOOD 

22.2 

6 

9159-9184 

GOOD 

10 

9185-9210 

GOOD 

20 

9211-9236 

GOOD 

40 

9237-9262 

GOOD 

60 

9263-9288 

GOOD 

27.8 

6 

8924-8949 

GOOD 

10 

8898-8923 

GOOD 

20 

8950-8975 

GOOD 

40 

8976-9001 

GOOD 

60 

9002-9028 

GOOD 

33.3 

6 

9289-9314 

GOOD 

10 

9315-9340 

GOOD 

XU 

9341-9366 

EXCELLENT 

9367-9392 

EXCELLENT 

uM 

9393-9419 

EXCELLENT 

37.0 

6 

9420-9445 

EXCELLENT 

10 

8798-8823 

GOOD 

20 

8824-8844 

GOOD 

40 

8845-8871 

GOOD 

60 

8872-8897 

GOOD 

51.8 

10 

9446-9471 

EXCELLENT 

20 

9472-9497 

EXCELLENT 

60 

9502-9527 

EXCELLENT 

120 

9530-9551 

EXCELLENT 
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